Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN M£X style file vl.4) 



Compressible Turbulence in Galaxy Clusters: Physics and 
Stochastic Particle Re-acceleration 

G. Brunetti, 1 A. Lazarian 2 

1 INAF - Istituto di Radioastronomia, via Gobetti 101, 1-40129 Bologna, Italy 

2 Department of Astronomy, University of Wisconsin at Madison, 5534 Sterling Hall, 475 North Charter Street, Madison, WI 53706, USA 

o 

Q ■ 5 February 2008 

H 

& ■ ABSTRACT 



(N 



OS 

in 

m 



We attempt to explain the non-thermal emission arising from galaxy clusters as a re- 
sult of the re-acceleration of electrons by compressible turbulence induced by cluster 
mergers. On the basis of the available observational facts we put forward a simplified 
model of turbulence in clusters of galaxies focusing our attention on the compress- 
ible motions. In our model intracluster medium (ICM) is represented by a high beta 
plasma in which turbulent motions are driven at large scales. The corresponding in- 
jection velocities are higher than the Alfven velocity. As a result, the turbulence is 
approximately isotropic up to the scale at which the turbulent velocity gets compa- 
rable with the Alfven velocity. These motions are most important for the energetic 
particle acceleration, but at the same time they are subjected to most of the plasma 
J*"-. , damping. Under the hypothesis that turbulence in the ICM is highly super- Alfvenic 

■ the magnetic field is passively advected and the field lines are bended on scales smaller 

^ (— | than that of the classical, unmagnetized, ion-ion mean free path. This affects ion dif- 

fusion and the strength of the effective viscosity. Under these conditions the bulk of 
turbulence in hot (5-10 keV temperature) galaxy clusters is likely to be dissipated 
£5 ■ at collisionless scales via resonant coupling with thermal and fast particles. We use 

r/} , collisionless physics to derive the amplitude of the different components of the energy 

of the compressible modes, and review and extend the treatment of plasma damping 
in the ICM. We calculate the acceleration of both protons and electrons taking into 
account both Transit Time Damping acceleration and non-resonant acceleration by 
large scale compressions. We find that relativistic electrons can be re-accelerated in 
the ICM up to energies of several GeV provided that the rms velocity of the com- 
pressible turbulent-eddies is (Vl/c s ) 2 w 0.15 — 0.3; c s is the sound speed in the ICM. 
We find that under typical conditions w 2-5 % of the energy flux of the cascading of 
compressible motions injected at large scales goes into the acceleration of fast particles 
and that this may explain the observed non-thermal emission from merging galaxy 
clusters. 

Key words: acceleration of particles - turbulence - radiation mechanisms: non- 
thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 



1 INTRODUCTION 

In the last years observations of galaxy clusters have shown that non-thermal components are mixed together with the thermal 
component of the intracluster medium (ICM) . A fraction of massive galaxy clusters hosts diffuse radio emission in the form of 
radio halos, Mpc-sized diffuse synchrotron radio sources at the cluster center, and radio relics, elongated diffuse synchrotron 
radio sources at the cluster periphery. This directly proves the presence of GeV relativistic electrons (and/or positrons) and 
of fiG magnetic fields diffused on Mpc scales through the cluster volume (e.g., Feretti 2005, for a review). A related issue 
is the discovery of non-thermal emission in the hard X-ray band detected in a few galaxy clusters (e.g., Fusco-Femiano et 
al. 2004; Rephaeli, Gruber, Arieli 2006). 
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The most spectacular example of non-thermal emission from galaxy clusters is that of giant radio halos. These are 
very extended (Mpc) synchrotron radio emissions, not connected with cluster radio sources, at the center of clusters, with 
a steep spectrum and a typical synchrotron luminosity in the range w 10 40 — 10 42 erg s _1 . A remarkable point is that the 
emitting particles have a life-time considerably shorther than that necessary to diffuse over the scales of these radio halos, 
and this poses a theoretical problem on their origin (e.g., Jaffe 1977). In principle if the content of cosmic ray hadrons in 
the ICM is sufficiently large, fast electrons and positrons may be continuously injected in the ICM via hadronic collisions 
between cosmic-rays and thermal protons (Dennison 1980; Blasi & Colafrancesco 1999), alternatively different forms of in 
siiw-stochastic acceleration and re-acceleration operating for a small fraction of the cluster life may provide a viable source for 
high energy emitting particles (Schlickeiser et al. 1987; Brunetti et al. 2001; Petrosian 2001). Future gamma ray observations 
(GLAST, Cerenkov telescopes) will constrain the content of cosmic-ray hadrons in galaxy clusters and provide an important 
tool to better understand the origin of the relativistic particles in the ICM (e.g. Reimer 2004 & Blasi et al. 2007 for recent 
reviews) . 

It is believed by several authors that the re-acceleration scenario may provide a promising picture to explain the bulk of 
present-day radio data (e.g., reviews by: Brunetti 2003,04; Petrosian 2003; Blasi 2004; Hwang 2004; Feretti 2005). This model 
essentially relies on the hypothesis that a fraction of the kinetic energy associated with cluster-cluster mergers is channelled 
into turbulence and re-acceleration of relativistic particles in the ICM. 

MHD turbulence is known to be an important agent for particle acceleration since Fermi (1949) first pointed this out. Second 
order Fermi acceleration by MHD turbulence was appealed for acceleration of particles in many astrophysical environments, 
e.g. Solar wind, Solar flares, ICM, gamma-ray bursts (see Schlickeiser & Miller 1998; Chandran 2003; Brunetti et al. 2004; 
Petrosian & Liu 2004; Cho & Lazarian 2006; Petrosian et al. 2006; Becker, Le & Dermer 2006; Dogiel et al. 2007). Naturally, 
properties of compressible MHD turbulence (see Shebalin, Matthaeus, & Montgomery 1983; Higdon 1984; Montgomery, Brown, 
& Matthaeus 1987; Shebalin & Montgomery 1988, Zank & Matthaeus 1992; Cho & Lazarian 2003 and references therein) are 
essential for understanding the acceleration mechanisms. 

Among the advances in understanding MHD turbulence, we would like to mention the Goldreich & Sridhar (1995; 
henceforth GS95) model of turbulence. GS95 dealt with incompressible MHD turbulence and showed that Alfven and pseudo- 
Alfven modes follow the scale-dependent anisotropy of Zy ~ f[ Z , where Zy is the size of the eddy along the local mean magnetic 
field and l± that of the eddy perpendicular to it. Lithwick & Goldreich (2001) conjectured that this scaling of incompressible 
modes is also true for Alfven modes and slow modes in the presence of compressibility. In Cho & Lazarian (2002; 2003, 
henceforth CL03) the rational for considering separately the evolution of slow, fast and Alfven mode cascades was justified. 
The numerical simulations in CL03 and Kowal & Lazarian (2006) verified that Alfven and slow mode velocity fluctuations 
are indeed consistent with the GS95 scaling, while fast modes exhibit isotropy in both gas-pressure (high /3 p i) and magnetic- 
pressure (low (3 p i) dominated plasmas. The former case is the most appropriate for clusters of galaxies that we deal in this 
paper. 

A single most important change in the paradigm that has become obvious recently is that if the turbulent-energy injection 
happens at large scales, the cascading Alfvenic mode is presented at sufficiently small scales by very elongated eddies. Thus the 
interactions of these mode with cosmic rays differs considerably from that of the isotropic eddies that earlier researchers dealt 
with. Under these conditions nearly isotropic fast modes were identified as the dominant agent for scattering and resonance 
acceleration (Yan & Lazarian 2002). As a consequence of this, our understanding of energetic particle-turbulence interactions 
via gyroresonance and the Transit-Time Damping (TTD) (Chandran 2000; Yan & Lazarian 2002, 2004; Farmer & Goldreich 
2004; Cassano & Brunetti 2005) as well as the non-resonance acceleration of cosmic rays by large scale compressible motions 
(see Ptuskin 1988, Chandran 2003, Cho & Lazarian 2006) has been altered. This calls for the corresponding advances in the 
treatment of cosmic-ray acceleration in the environment of clusters of galaxies (see e.g. Brunetti 2006 & Lazarian 2006a) . 



2 OUTLINE OF THE PAPER 

In this paper we proceed in three main steps : 

I) As a first point we discuss the problem of turbulence in the ICM and work up a simplified but physical picture of the 
properties and relevant scales of turbulence in galaxy clusters. As we discuss below (see §3.1) the plasma in clusters of galaxies 
is expected to be both magnetized and turbulent. Its Reynolds numbers are expected to vary as the magnetic field grow 
(§3.2), but they are expected to be high enough to allow turbulence to be excited. Finally turbulence in hot galaxy clusters 
is expected to be dissipated via collisionless dampings and this makes the particle acceleration process a natural consequence 
(§3.4). This part of the paper is mainly designed to provide a reference picture for observers and a viable astrophysical starting 
point for theoretical developments. 

II) As a second point we discuss the physics of compressible motions in the collisionless regime. This part of the paper 
is a necessary extension of previous seminal studies of collisionless turbulence and is aimed at the presentation of necessary 
general equations to use in the paper. In particular to characterize the plasma-cosmic rays interactions we characterize the 
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compressible motions using dielectric tensor (§4.1), give the expression for the energy spectrum of the fast modes in §4.2 and 
describe the TTD damping in intracluster plasma in §4.3; complex expressions and calculations are reported in Appendices 



III) Finally we discuss the issue of stochastic particle acceleration in galaxy clusters by compressible turbulence. The 
resonant TTD acceleration is discussed in §5.1, while the effect of non-resonant acceleration is discussed in §5.2. In §6 we 
discuss the results in the framework of the particle re-acceleration model in galaxy clusters: in §6.1 we briefly review the 
basic physics of cosmic rays in galaxy clusters, and in §6.2 we present detailed calculations on particle re-acceleration in the 
ICM. Here we claim that compressible turbulence may drive efficient particle acceleration in the ICM. This is an extension of 
recent studies on the argument and provides a view of the process of particle re-acceleration in the ICM which is additional 
(or alternative) to that of Alfvcnic acceleration. 

In §7 we discuss the most relevant findings and simplifications, and in §8 we provide a short summary. 

3 TURBULENCE IN THE ICM 

3.1 Injection of turbulence in the ICM 

Cluster mergers and accretion of matter at the virial radius may induce large-scale motions with Vl ~ 1000 km s _1 in massive 
clusters. Numerical simulations suggest that turbulent motions may store an appreciable fraction, 5-30%, of the thermal energy 
of the ICM (e.g., Sunyaev, Bryan & Norman 2003; Dolag et al. 2005; Vazza et al. 2006). Simulations of merging clusters provide 
an insight into the gas dynamics during a merger event (e.g., Roettiger, Burns, Loken 1996; Roettiger, Loken, Burns 1997; 
Ricker & Sarazin 2001): sub-clusters generate laminar bulk flows through the sweeped volume of the main clusters which 
inject turbulence via e.g. Kelvin-Helmholtz instabilities at the interface of the bulk flows and the primary cluster gas. The 
largest turbulent eddies decay into smaller and turbulent velocity fields developing a turbulent cascade. 

A simple, but well motivated by physical arguments, semi-analytical approach allows to follow the cosmological injection 
of merger-turbulence. Calculations from Cassano & Brunetti (2005) suggest that turbulence in the ICM is transient being 
mostly injected during the most massive mergers. However, since more frequent minor mergers may also contribute to the 
injection of such turbulence, some minimum level of turbulence should be rather ubiquitous in the ICM. In these calculations 
turbulence is assumed to be injected in the cluster volume swept by the sub-clusters, which is bound by the effect of the ram 
pressure stripping, and the turbulent energy is calculated as a fraction of the PdV work done by the sub-clusters infalling 
onto the main cluster. Essentially merger-driven turbulence is powered by the gravitational potential well and thus the energy 
of this turbulence should approximatively scale with the cluster thermal energy (Cassano & Brunetti 2005) . Support to this 
scaling comes from a recent analysis of a sample of galaxy clusters from cosmological numerical simulations (Vazza et al. 2006). 

Turbulence is an important ingredient in the physics of the ICM as it is necessary to understand the amplification of 
magnetic fields in clusters (Dolag et al. 2002; Schekochihin et al. 2005; Subramanian et al. 2006), an issue closely related 
to the non-thermal emission from clusters but that we will not address in this paper. Turbulence might provide a source of 
heating to balance the cooling of cluster cores (Fujita, Matsumoto & Wada 2004), and the knowledge of the basic aspects 
of turbulence in galaxy clusters is also crucial to model the transport of heat and metals in the ICM (Narayan & Medvedev 
2001; Cho et al., 2003; Voigt & Fabian 2004; Lazarian 2006b). 

In spite of obvious observational challenges, indications of some level (at least 10-20% of the thermal energy) of turbulence 
in the ICM comes from gas-pressure maps in the X-rays (Schuecker et al. 2004) , and also from the lack of resonant scattering 
from X-ray spectra (Churazov et al. 2004; Gastaldello & Molendi 2004). 

Interestingly enough, also upper limits to the turbulent -energy content in the ICM were obtained in a few nearby galaxy 
clusters from kinematical arguments related to the properties of Ha and X-ray filaments (e.g. Fabian et al. 2003; Crawford 
et al. 2005; Sun et al. 2006). Assuming that turbulence is driven at hundred-kpc scales the above upper limits actually can 
be used to place upper limits on the intensity of strong turbulence in the ICM (supersonic or trans-sonic turbulence). 

3.2 Reynolds Number and developing of turbulence in the ICM 

In this Section we discuss the important issue of the Reynolds number of the fluid in the ICM, and derive its value by assuming 
a simple, but physically motivated, scenario. 

A fluid becomes turbulent when the rate of viscous dissipation at the injection scale, L , is much smaller than the energy 
transfer rate, i.e. when the Reynolds number is Re — VlL /vk 1, where Vl is the injection velocity and vk is the kinetic 
fluid viscosity. The main source of uncertainty here comes from our ignorance of vk in the ICM. 

If the ICM were not magnetized vk ~ lmf P Vi/3, were Vi is the velocity of thermal ions, and l m f P is the ion-ion mean free 
path in case of pure Coulomb interactions (e.g., Braginskii 1965): 



A-C . 




(1) 
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where In A is the Coulomb logarithm. 

Thus the corresponding Reynolds number would be: 

Re gnf — \( L ° \( Uth \( T Y 5/2 (—\ (2) 
Vl000km/s/ V300kpc/ U0- s cm- 3 / V8kcV/ V 40 / ( > 

which is formally just sufficient for initiating the developing of turbulence. 

However, in the presence of (even a small) stationary magnetic field the Reynolds number for motions in the direction 
perpendicular to the magnetic field gets extremely high essentially because the perpendicular mean free path of particles is 
limited to the Larmor gyroradius-scale (e.g., Braginskii 1965). Potentially, diffusion along the wandering turbulent-magnetic 
field lines could significantly increase the particle diffusivity and the plasma viscosity. For instance, estimates in Narayan & 
Medvedev (2001) suggest that electron diffusivity in a turbulent medium can be of the order of 1/5 of the classical Spitzer 
value for unmagnetized medium, provided that the injection velocity, Vl , is equal to the Alven velocity. 

More general calculations (Lazarian 2006a, and ref. therein) show that things could be more complicated and that the effective 
viscosity depends on the super- or sub- Alfvenic nature of the turbulence*. Turbulence in the ICM is super- Alfvenic, i.e. 
turbulence with the injection velocity larger than the Alfven one. In this case the turbulent hydrodynamic motions can easily 
bend the magnetic field lines. The trajectory of the particle that follows such a field line gets diffusive even in the absence of 
collisions. The effective mean free path of a particle is determined by the scale at which magnetic tension can withstand the 
hydrodynamic forces, i.e. the scale at which the turbulent velocity, Vi A , gets equal to the Alfven one, va = B/^/4-wp, where 
B ~ (Bo + Brms) 1 ^ 2 , (Lazarian 2006b). This scale, at which turbulence becomes MHD, is t Ia ~ L (^) 3 and assuming 
typical conditions in Mpc regions at the center of massive (and hot) galaxy clusters where radio halos are found, it is : 



which is < l m fp- This implies the important point that, even for motions along the magnetic field, the Reynolds number in a 
turbulent ICM is larger than that estimated with the classical formula (Eq.2). Actually one finds Re > few times 10 3 which 
ensures that the ICM gets turbulent. 

The main uncertainty in the evaluation of the Reynolds number comes from the value of the effective mean free path of 
particles. Eq.(3) accounts for the effect of the turbulent magnetic field, however additional mechanisms may affect the value of 
the particle mean free path in the ICM, for example plasma instabilities. Plasma instabilities could be at work in the ICM, e.g. 
turbulent compressions themselves may drive instabilities. These instabilities in the ICM may induce scatterings of thermal 
ions which reduce the effective mean free path and further increase the value of the Reynolds number (e.g., Schekochihin & 
Cowley 2006; Lazarian & Beresnyak 2006). In what follows to be conservative and with the aim to simplify the overall picture, 
we disregard this effect, so that our estimates of the Reynolds number in the ICM would be considered as a lower limit. 

3.3 Turbulent Modes 

3. 3. 1 Basic properties of the turbulent modes 

Turbulence discussed in the previous Sections is a complex mixture of several turbulent modes. The ICM is a compressible 
high-beta plasma. At large scales, where magnetic fields are not dynamically important (Vl >> va), the turbulence is 
essentially in the hydro- regime, and we shall assume that turbulence in the ICM is done by solenoidal and compressible 
(essentially sound waves) motions. At smaller scales, in the MHD regime, it is Vi < va and three types of modes should 
exist in a compressible magnetized plasma: Alfven, slow and fast modes. Slow and fast modes may be roughly thought as 
the MHD counterpart of the compressible modes, while Alfven modes may be thought as the MHD counterpart of solenoidal 
Kolmogorov eddies (a more extended discussion can be found in Cho, Lazarian & Vishniac 2002, and ref. therein). Sound 
modes at large scales have propagation properties similar to that of the fast modes in the MHD- regime. For this reason in 
this paper we shall use the properties of these modes for describing compressible turbulence, i.e. hydro- modes (sound waves) 
at large scales and fast modes themselves at small scales. Fast modes are compressive waves which propagate across or at an 
angle to the local magnetic field. The fast mode branch in a plasma extends from low frequencies up to the electron cyclotron 
frequency. At frequencies below the ion cyclotron frequency, and in the weak damping limit, the dispersion relation of these 
modes is given by lo = V p hk, where the phase velocity is given by (e.g., Krall & Trivelpiece 1973) : 



>=|(/W2 + l))l + ^-4(^) 



Vph s ° 4 i k j (d+«i) 2 f 1 i ~ ' \r "w a+/w2) 



* The super- and sub-Alfvenic are determined in terms of the total magnetic field, 
t In deriving Eq.(3) we use the hydro- scaling V; oc i 1 / 3 
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and where the parameter beta of the plasma is defined by /3 p i — 2c 2 /v\. 

Alfven modes propagate along or at an angle to the local magnetic field. The Alfven branch extends from low frequencies 
up to the ion cyclotron frequency. In this frequency range the Alfvenic dispersion relation is given by u> = UA|fc|||, where 
va = B /^/4irn th mi is the Alfven velocity. 

Alfven and fast modes differ also for the direction of the displacement vectors: the displacement of Alfven modes is 
always perpendicular to B , while that of fast modes makes an angle to the local magnetic field and in the case fi p i — > it 
is perpendicular to B , while for (3 p i — > oo it becomes radial (along k); a detailed discussion on the decomposition of MHD 
modes can be found in Cho & Lazarian (2002), (2003) and in Kowal & Lazarian (2006). 

Slow modes has "-" before the square root in Eq.(4) and for /3 p i » 1 they have the dispersion relation of Alfven modes. 
We will not include slow modes in our calculations in the in the particle acceleration process by large scale modes (Sect. 5) 
as they have a phase velocity in the ICM much smaller than that of the fast modes and thus are less important. 

At MHD-scales Alfven and slow modes might be of some relevance in discussing the particle acceleration process either 
because they can accelerate particles, or because in principle they may provide a particle pitch-angle scattering process^ 
which is required by acceleration processes driven by other modes (§ 4,5). 

3.3.2 Coupling between turbulent modes in the ICM 

Although the complex dynamics of galaxy clusters and the relatively large value of the Reynolds number of the ICM are likely 
to make the ICM itself a turbulent medium, it is somewhat difficult to have a clear idea of the relative importance of the 
different turbulent modes in the ICM. Indeed this depends on the nature of the turbulent forcing and on the mode coupling 
between different modes in the ICM. 

We shall assume that a sizeable part of turbulence at large scales (namely at scales where the magnetic tension does 
not affect the turbulent motions) is in the form of compressible motions. This is reasonable as these modes are expected to 
be easily generated in a high beta medium even in the unfavourable case of solenoidal turbulent forcing. This is proved by 
closure calculations carried out in the case of f3 p i >> 1. Indeed when motions are hydro- in nature the coupling between 
solenoidal and compressible motions is efficient and the excitation of compressible modes by the solenoidal modes is driven by 
the incompressible pressure arising from non-linear interaction between solenoidal modes themselves (Bertoglio et al. 2001). 
These studies have shown that the fraction of energy in the form of compressible modes is found to scale with oc M 2 x Re for 
M 2 x Re < 10 (M s < 1 is the turbulent Mach number), while for M 2 x Re > 10 the scaling is expected to flatten (Bertoglio 
et al. 2001; Zank & Matthaeus 1993). Obviously a solenoidal turbulent forcing, which limits the energy of compressible modes 
to be smaller than that of solenoidal modes (even in the super-Alfvenic case) , is probably not appropriate for galaxy clusters 
where turbulence is likely to be excited by compressible forcings, and this might result in a larger ratio between compressive 
and solenoidal modes (at least for super-Alfvenic motions). 

Situation may be radically different at smaller scales where the magnetic field tension affects turbulent motions, i.e. in the 
MHD- regime, I < I a- In this case, MHD numerical simulations have shown that a solenoidal turbulent forcing gets the ratio 
between the amplitude of Alfven and fast modes in the form (Cho & Lazarian 2003) : 

{SV)l (3V)sVA ... 
{5Vf s ~ cj + v\ w 

which essentially means that coupling between these two modes may be important only at I ~ I a (in the MHD- regime it 
should be [SV) S < va) since the drain of energy from Alfvenic cascade is marginal when the amplitudes of perturbations 
become weaker. Most importantly in galaxy clusters it is c 2 >> v\ and thus the ratio between the amplitude of Alfven and 
fast modes at scales I < Ia is expected to be small, (SV) 2 /(SV) 2 < (va/c s ) 2 ~ 10~ 2 (this for solenoidal forcing at I « Ia)- 

3.4 Dissipation of turbulence in the ICM 

3.4-1 Collisional regime & viscous dissipation 

Viscosity is important in the dissipation of turbulent eddies in the collisional regime. In this regime the cascade of hydro- 
motions is maintained down to a scale Idiss ~ L (Re)~4 at which the viscous dissipation rate equals the wave energy transfer 
rate. The damping rate of hydro- motions at scale I due to the viscosity is : 

n ~ v f (6) 

t The Alfvenic mode as well as the slow mode gets anisotropic for scales less than I a, which makes the scattering inefficient, however 
(Chandran 2000, Yan & Lazarian 2002). 



© 0000 RAS, MNRAS 000, 000-000 



6 G. Brunetti, A. Lazarian 



here vk is a reference value of the kinetic viscosity which gives the main uncertainty in the calculations. 

As a simplified and conservative approach we can assume that B D is initially ordered and that the first super-Alfvenic 
turbulent eddies, injected at large L » l m fp scales, initiate a cascading and thus that the bending of the field lines follows 
this cascading. Turbulent motions along B experience the strongest viscous dissipation which can be grossly estimated by 
using the classical formulation of (unmagnetised) viscosity. By taking physical conditions appropriate for the central Mpc of 
hot galaxy clusters, the dissipation scale of these parallel motions reads : 

/ ~; ( Vi \ l ( L °\*~l ( Vl Y~ 4 ft L ° w v 8keV i InAA* 

diss ~ mfp \W L J \lmSv) ~ Vl03km/ S J ^300kpc n iO-3 K T > { 40 } ) K<) 

while turbulent motions transverse to B experience a much smaller viscosity and shall cascade at scales << l m f P - The 
cascading of these transverse (quasi-perpendicular) motions at a given scale takes a time of the order of the the bending time 
scale of the magnetic field on the same scale and these motions become the responsible for the bending of the field lines on 
scales << l m fp, potentially down to scales « I a- We note that indeed recent Bayesian analysis of RM show that magnetic 
fields in galaxy clusters could be tangled at least on scales « kpc (Vogt & Ensslin 2005), smaller scales being inaccessible to 
observations, thus suggesting that the bending of the field lines happens on scales << l m fp- 

As discussed in § 3.2 the bending of the field lines on scales < l m f P reduces the effective particle mean free path yielding a 
reduction of the viscosity. Viscosity indeed depends on the flux of the momentum which is transported by particles and this 
is determined by the diffusion of the particles that carry this momentum from the layers moving with different velocities. By 
limiting this diffusion the turbulent-bending of the field lines decreases the viscosity and thus the dissipation of turbulence 
itself. 

Thus the turbulent eddies which cascade afterwards evolve in a very tangled magnetic field and experience an effective viscosity 
which we shall adopt in the form vk « 1/3v»Za, and the effective dissipation scale, l b diss , we would encounter in the case of 
super-Alfvenic turbulent ICM becomes : 

^ lss ~ ' dlsa (l^) ~ 45 lmfp (f0«km/s) ( 300k^c" } (^) (lO-wO (skcv) ( lfT ) (8) 

Also in this case the effect driven by plasma instabilities in the ICM may affect our estimates. In particular, the scattering 
of thermal ions induced by these instabilities may additionally decrease the effective viscosity in the ICM, and this might 
reinforce our conclusions that, even assuming collisional physics, the bulk of compressible turbulent motions in the ICM is 
expected to be dissipated only at small scales, < kpc. 



3.4-2 Collisionless regime 

The viscous damping is not important in the collisionless regime, i.e. when the scales of interest are smaller than the particle's 
mean free path or when the time-scales of interest are shorter than the particle's collision time. When the diffusive-trajectory 
of particles is not driven only by collisions (as indeed in the super-Alfvenic turbulent-magnetized case, § 3.2, 3.4.1) the most 
appropriate way to define the collisionless regime is in terms of collision frequency, and we shall use collisionless physics for 
the turbulent modes when the frequency of these modes is larger than the ion-ion collision frequency uu (e.g., Eilek 1979) : 

4 r e 4 n th ln(A) 

m * -^ m y\ kBT yn (9) 

Magnetosonic modes dissipate energy in the collisionless regime in accelerating charged particles especially via Transit-Time- 
Damping (e.g., Schlickeiser & Miller 1998) which is particularly severe in high beta-plasma conditions like those in the ICM. In 
terms of scales, from Eq.(9) and uj — V p hk, the collisionless regime for magnetosonic waves in the ICM starts approximatively 
at the scale of the ion mean free path lc ~ lmfp (Eq. 1). Thus from a general point of view, in order to understand the way 
compressible modes dissipate in the ICM it is necessary to compare the viscous dissipation scale, l b diss , with the collisionless 
scale lc'- if l b diss < lc the cascading process of these turbulent modes would reach collisionless scales before being significantly 
affected by viscosity and energy will be dissipated via collisionless dampings, while in the opposite case turbulence will be 
dissipated by viscosity. 

From Eqs.(7-8) we immediately have that the bulk of compressive turbulence in the hot ICM is likely to be dissipated via 
collisionless dampings. Indeed in hot (and massive) galaxy clusters it is found that viscosity is not efficient enough to dissipate 
the turbulent motions, unless the large-scale velocity of these motions is relatively small, Vl < 300 km/s, namely in case 
of very low turbulence. At the same time, however, an efficient dissipation of turbulent motions may happen in strongly 
magnetized (B > 5/iG) , lower temperature and high density regions which are conditions appropriate at the center of clusters 
with cooling flows (cool cores). Here viscosity may potentially become an important source of dissipation of the turbulent 
eddies. 
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GC 


CC 


galactic halo 


HIM 


WIM 


Sun 


T(K) 


10 8 


3 x 10 7 


2 X 10 6 


10 6 


8 x 10 3 


10 7 


c s (km/s) 


1650 


900 


130 


90 


8 


360 


n th (cm~ 3 ) 


10~ 3 


5 x 10~ 2 


10~ 3 


4 x 10~ 3 


0.1 


10 10 


Imfp ( cm ) 


5 x 10 22 


10 20 


4 X 10 19 


2 x 10 18 


6 x 10 12 


10 8 


Lo(pc) 


1 - 5 x 10 5 


1 - 5 x 10 5 


100 


100 


50 


3 x 10" 10 


B(juG) 


1 


10 


5 


2 


5 


10 8 


ciK 


500 


100 


0.3 


3.5 


0.1 


0.03 


damping 


collisionless* 


collisionless ? 


collisionless 


collisional 


collisional 


collisionless** 



Table 1. The reference parameters of astrophysical plasma and relevant damping. The dominant damping mechanism for turbulence 
is given in the last line. GC= hot galaxy clusters (* Vl > 300 km/s), CC= cool-cluster cores, HIM= hot ionized ISM, WIM= warm 
ionized ISM, SUN= Solar flare plasma (** Alfvcnic turbulcnt-Mach number Ma > 0.3 is assumed). 

It should be mentioned that plasma instabilities might complicate the picture. On one hand, their straightforward effect is 
to decrease the effective viscosity in the ICM, however, on the other hand they introduce a new relevant scattering frequency 
of ions which could be larger than the ion-ion scattering frequency (Eq. 9) and the net result might be that the collisionless 
regime gets into play at smaller scales. As in the previous Sections we discard this possible effect which would deserve detailed 
investigation. 

The nature (collisional or collisionless) of the turbulent dissipation in astrophysical plasma is a crucial point. In Tab. 1 we 
report the case for several astrophysical situations undertaking different physical conditions, processes and scales of interest. 
A collisionless dissipation of compressive turbulence is believed to be eventually operating in a few other astrophysical regions 
such as in solar flare plasma and in the Galactic Halo. It is important to note here that stochastic particle acceleration is 
indeed suggested to power the hard X-ray flares observed in the Sun (e.g., Miller, La Rosa & Moore 1996; Petrosian, Yan & 
Lazarian 2006). We note that the beta of plasma in these filaments is extremely small and thus even in the case of strong 
turbulence the collisionless dissipation of compressive modes should happen at scales, I < I a, at which turbulent motions are 
MHD in nature. On the other hand turbulence in the ICM is super-Alfvenic (essentially due to the high beta of plasma) and 
the collisionless regime in the hot ICM starts at scale I > Ia were compressive motions are still hydro- in nature. 

3.5 Conclusion I: Turbulent Scenario in the ICM 

Given the above discussions it is possible to set up a simplified and operative scenario of turbulence in the ICM to adopt in 
this paper. 

Within a simplified picture of turbulence that we consider here, super-Alfvenic turbulence is made by a mix of magnetosonic 
modes (essentially similar to sound modes) and incompressible-Kolmogorov turbulent eddies (which roughly correspond to 
the Alfven modes in the MHD regime). 

We shall assume that turbulence is injected at large scales L ~ 300 — 500 kpc most likely by a complex mixture of compressive 
and solenoidal forcing. The typical velocity of the turbulent eddies at the injection scale is expected to be around Vl ~ 
500 — 1000 km/s which makes turbulence sub-sonic, with M s — Vl/c s w 0.3 — 0.8, but strongly super-Alfvenic, with 
Ma = Vl/va > 10. Turbulent motions at large scales are thus essentially hydrodynamics and the cascading of compressive 
(magnetosonic) modes may couple with that of solenoidal motions (Kolmogorov eddies) . 

Assuming typical conditions in hot (and massive) galaxy clusters we find that even in the unmagnetized case viscosity would 
still allow hydro- motions to cascade down to scales of the order of « Imfp- In the magnetized case viscosity is believed to be 
partially suppressed. In addition when turbulence is super-Alfvenic hydro- motions can easily bend the magnetic field lines 
affecting the effective mean free path of ions which happens to become limited approximatively to the MHD scale, Ia§ The 
value of the effective viscosity, even for motions along the magnetic field lines, is thus expected to be considerably reduced 
with respect to the classical unmagnetized value and one may adopt a reasonable value of the Reynolds number Re > 10 3 . 
The important consequence of this picture is that both solenoidal and compressive modes in hot galaxy clusters would not be 
strongly affected by viscosity at large scales and an inertial range is established, provided that the velocity of the eddies at 
large scales exceeds « 300 km/s. We shall assume that a sizeable part of the large scale turbulence is done by magnetosonic 
(essentially sound) modes. At collisionless scales, / < 10 — 50 kpc, these modes are affected by strong collisionless dampings 

§ This provided that turbulent eddies may reach the MHD scale without being dissipated (§ 5.1.3, § 5) 
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with both thermal and relativistic particles (§ 4.3) and thus they are expected to be the modes which dominate the particle 
acceleration process. Our claim about the existence of this well developed turbulent cascade which establishes an inertial 
range from large scales to the collisionless scales would be even reinforced when the possible effect of plasma instabilities is 
considered. 

Although in this paper we focus on the particle acceleration by hydro- magnetosonic modes, it is worth mentioning that the 
mode composition at smaller scales, I « I a, in the ICM should becomes much complex. We shall assume that Alfven modes 
are present at these MHD scales in the ICM since in principle the cascading of solenoidal modes might reach very small 
scales, due to the lack of large-scale collisionless dampings for these modes, and also because several mechanisms can convert 
a fraction of the energy flux of large-scale turbulent cascade in the injection of Alfven modes at smaller scales (e.g. Kato 
1968; Eilek & Henriksen 1984; Lazarian & Bercsnyak 2006). At scales I < I a, the coupling between Alfven and compressible 
modes gets changed and only slow modes are cascaded by Alfvenic modes (GS95, Lithwick & Goldreich 2001, Cho & Lazarian 
2002), while the cascading of fast modes is not particularly sensitive to the presence of the other modes. Given that, and 
since magnetosonic modes are expected to be damped at scales larger than (or similar to) I a (§ 4,5), the spectrum of the 
ICM-turbulence at I << I a is expected to be populated only by Alfven and slow modes. These modes however would get 
anisotropic at these scales (unless injected at these scales by some mechanism) and this should reduce their contribution to 
the scattering and acceleration of fast particles via gyro-resonance. 



4 COMPRESSIBLE TURBULENCE IN THE COLLISIONLESS REGIME 

Compressible turbulence in galaxy clusters is thus made of large scale hydro-motions with frequencies essentially infinitely 
small with respect to S~li/f3 p i (fi; being the Larmor frequency of ions). The basic physics of these low-frequency compressible 
modes in the collisionless regime can be derived by mean of the quasi-linear theory and has been investigated in several 
seminal papers (e.g., Melrose 1968; Barnes 1968; Baldwin, Bernstein & Weenink 1969; Barnes & Scargle 1973, hereafter 
BS73)^. This Section extends previous studies as we derive specific and operative expressions for the physical properties of 
these modes which are of interest for the present paper (e.g., energy decomposition of the mode, TTD-damping rate) and 
discuss their dependence on the mode-propagation angle. We focus on the case of long-wavelength modes in a magnetized 
plasma dominated by thermal particles as it should be the case of the ICM. Here we report the main formulae, while details 
and derivation of the main equations are given in the Appendices. 



4.1 Geometry of the Mode and Dielectric Tensor 

We define the turbulent fluctuations associated with the electric and magnetic field as : 

E — TZ (Efc exp{i(k • r — cut)}) (10) 
and 

B = ft(B fc cxp{i(k-r-a;t)}) (11) 

where HQ stands for the Real part. In the collisionless regime it is usual to start with fixing the geometry of the mode 
propagation and of the electric field fluctuations. Without loss of generality we may chose the particular system where the 
^/-component of the wavevector of the modes vanishes, i.e. : 

k=(* ± ,0,fc|l) (12) 

For this choice the amplitude of the electric field (and spatial Fourier transform of the electric field of the mode) is given by 
(e.g., BS73) : 

E fc = (0,B±,S||) (13) 
The amplitude of the magnetic field of the mode comes from the Faraday low, B t = x Et : 

B k = -(-k^E ± ,-k±E h k±E ± ) (14) 



c 

As a starting point we assume the presence of several, a, species of particles with particle momentum given by 



Pt. = (pxcos0,p x sin0,p||) = m Q 7 (y± cos<p,v ± sinewy) (15) 



^ For hydromagnetic waves with frequency of the order of Qi//3 p i see Footc & Kulsrud (1979). 
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and indicate with f a the normalized particle distribution in the momentum space of species a (f a (p) = N a f a (p)). 
The properties of a wave propagating in a magnetized plasma in the collisionless regime depend on the dielectric tensor of the 
plasma. In the general case, the dielectric tensor of the magnetized plasma is given by (e.g., Melrose 1968; see also Tsytovich 
1977 for the unmagnetized case) : 



Kij — Sij + 2n m G 



ni / _d_ 

■ ' l V ± dp\\ 



d 



{ : 



to — nil a - 



fell f ii 



k\\v\\ d 



(16) 



'± dp± dp 

where a» p , Q = \J 4-KN a e%/m a is the plasma frequency for the species a, hi — ^y§^ is the unit vector along the magnetic 
field, Q a = (e a B /m a c)/~{ is the Larmor frequency of particles a, 



V 



-i v -^J n (z)J' n {z) vl(JUz) 
V -^Jl{z) 



— iv±V\\ J n {z)J' n (z) 

iv±_v\\J n (z)J' n (z) v\J*(z) J 



(17) 



and z a = k±p±/m a Q" (f2" = Q'y is the classical Larmor frequency) is an adimensional parameter which scales with the ratio 
between the frequency of the mode and the particle Larmor frequency, and with the ratio between the particle velocity and 
the phase velocity of the mode, z a w (Lu/fl")(v a /V p h). We notice that magnetosonic modes with long wavelength, I >>pc, 
always have z a << 1. In this case a more suitable expression for the dielectric tensor can be obtained by expanding the 
Bessel-functions in Eqs. 16-17 in the limit z a << 1 (Appendix A). 



4.2 Energy of the Mode 

The energy of a mode in a magnetized plasma is done by the sum of the energy associated with the fluctuations of the electric 
and magnetic fields, We and Wb, and by the energy contributed by particles to the modes, Wp. The total energy is then : 

W = W B + W E + W P (18) 

In the collisional regime (and adiabatic equation of state) Wp is given by the contributions from the kinetic energy, Wk , and 
from a potential energy, W<s>, associated with pressure fluctuations, and a simple equipartition condition exists (e.g., Denisse 
& Delcroix 1963; Melrose 1968), namely: 

W B + Wt ^We + Wk (19) 

In the collisionless regime the medium is described in terms of the dielectric tensor and it is not possible to define Wp in a 
meaningful way. Thus one has to use Eq.(18) as a definition for Wp, with the total energy, W, defined independently. The 
total energy of the mode is given by (e.g., Barnes 1968; see also Melrose 1968 & Tsytovich 1972 for equivalent expressions): 

W(k,u) = ^[BUBki + E k 'i-^(uK^E k ^ (20) 

where stands for the Hermitian part of the dielectric tensor. In this case the first term in Eq.(20) accounts for magnetic 
field fluctuations, while (from Eq. 18) the term 

i^(^£H)^) = ^r + E(^))« (2D 

accounts for the contribution to the mode energy from the electric field fluctuations and from particles (Barnes 1968). 

In the quasi-linear regime the energy of the magnetic field fluctuations is related to that of electric field fluctuations by (e.g., 

Melrose 1968): 

= (c V |\ k E k 



WB = [yz) V-h-mn^ (22) 

which can be taken Wb — (c/V p h) 2 WE since under the physical conditions of interest for this paper magnetosonic modes have 
E±/E\\ << 1 (Appendix B). Thus combining Eq.(18) with Eqs. (20-22) it is easy to get the ratio between the total energy in 
the mode and that associated with the different components, Wb, We, and Wp. 

Thermal particles in the ICM should provide the dominant contribution to the total energy of turbulent modes. Thus we 
make the approximation that the dielectric tensor of the ICM is described by that of an electron-proton magnetized plasma 
in thermal equilibrium. Assuming a Maxwellian distribution for the thermal electrons and protons in the ICM : 
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Figure 1. Panel a) The expression S{j3 p i,9) is given as a function of 8. The behavior around 8 = tt/2 is highlighted in the right 
panel. kgT =8.6 keV is assumed. Calculations are reported for : c 2 s /v 2 A =ft p i/2= 100 (solid line), 10 (dotted line), 1 (dashed line), 0.5 
(long-dashed line), and 0.1 (dash-dotted line). 



up) = aw = (27m J; r) a /2 exp { ~ ^t? } (23) 

in Appendix B, from Eqs. (16-17) and Eq.(20) we show that the total energy of a fast mode is : 

™-^**[(£AK*)V^)*i(*)'(lM]> 

where the function S(f3 p i,9; {frei(p),T}) (Fig. 1) accounts for the terms of the dielectric tensor in the form : 



/ 



d Pndp±P ?l^ (25) 



which all come from the collisionless resonance between particles and modes with n — in Eq.(16) (see § 5.1 & Appendix 
B). Provided f a is an even function of py, only particles with velocity larger than the phase velocity of the mode can 
contribute to 5, since they should satisfy the condition fc|| «|| ~ w = V p hk. The velocity of the selected resonant particles scales 
as v ~ V p h{/3pi)/ cos(#), thus with increasing the angle between k and B , 9, particles with increasing velocities may contribute 
to this term. Formally particles with v — > oo contribute to S for 9 — > n/2, and this gets 5 — > in this limit (Fig. 1). Two 
wave-like behavior of S can be recognized in Fig. 1: the first one, for 6 < 1, marks the contribution from protons, while the 
second one, for larger 8, marks that from electrons, which are faster than protons. The resonance condition, kuv« ~ u>, also 
drives the shift of these wave-like behavior toward smaller values of 6 with decreasing f3 p i : when f3 p i decreases the resonance 
between the mode and a fixed portion of the particle distribution comes up at smaller values of 6. 

In Fig. 2 we report the ratio between magnetic and total energy of a mode propagating at an angle 6 as a function of f3 p i ; 
this ratio is independent of the wavenumber k of the mode. Two wave-like behaviors (due to the contribution from 5-terms) 
are visible: the first one, for 8 < 1, marks the contribution from protons, and the second one, for larger 9, from electrons. For 
small values of f3 p i it is Wb ~ W/2 and thus the quantity ft p i \Bk\ 2 /2W oc (3 p i, on the other hand, for large values of f3 p i it is 
Wb oc W/f3 p i and f3 p i\B k \ 2 /2W becomes independent of f3 p i. 

Finally, in Fig. 3 we report the ratio between particle energy and magnetic energy in the mode for different values of j3 p i (see 
caption). For (3 p i « 1, Wp reaches equipartition with Wb similarly to the case of collisional and low f3 p i plasmas. 
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0(rad) 0(rad) 



Figure 2. The ratio between the magnetic energy density of the mode and the total energy density of the mode is given as a function 
of 6 (for a better comparison the quantity is multiplied by 16irf3 p i/2). The behavior around 8 = tt/2 is highlighted in the right panel. 
kgT =8.6 keV is assumed. Calculations are reported for: c^,/v 2 A =f) p i/2= 100 (solid line), 10 (dotted line), 3 (dashed line), 1 (long— dashed 
line), 0.5 (short-dashed-dotted line), and 0.1 (long-dashed-dotted line). 



4.3 Turbulence Damping: TTD— resonance (n=0) 

A compressible mode in the collisionless regime experiences strong collisionless damping with thermal and relativistic particles 
and gets modified. In this Section we report relevant formulae for the collisionless damping rate via TTD-resonance of 
magnetosonic waves which will be used in the present paper (§ 5). 

The damping coefficient of the mode can be obtained by the standard formula for the linear growth rate of the mode in the 
quasi-linear theory (e.g., BS73)H; 

r = -<i^Lo^ (26) 

where A"?- stands for the anti-Hermitian part of the dielectric tensor, and u) r is the real part of ui. The general formula for 
the collisionless damping rate (n = 0, ±1, ..) is (Appendix C, and BS73): 



tv! a\ 71 fc ll 2 F , F , 2,,,af / w \ df a (p) df a (p)\ x ( P\\ , nn 0}a ~uj r j 



where 

2 



=2 



',J^(z a )Ej_ + — J n (z a )En 

p± 



(27) 



(28) 



In this paper we focus on the case n = (Transit Time Damping, discussed in § 5) which is the most important collisionless 
resonance between magnetosonic waves and particles in the ICM. In the case of long-wavelength magnetosonic waves the 
damping rate due to TTD-resonance with thermal electrons and protons with number density N e / p , is given by (Appendix 
C): 

r (k f>)- F |i3fc|2 -hU Vph M V ' h ( k \( k± Y ( m °"> kBT ) 1/ \ r Pxn l m */J& ( fc / fc ") 2 ( o Q) 

v ^°)-\l^w(kj) y ~ w^\M)Vk) i-c-^y N ^ exp {~ 2k B r i- ( ^ )2 r (29) 

where H(x) is the Heaviside step function (1 for x > 0, and otherwise), and the ratio |Bfc| jW is given by Eq.(24). 
Actually for a fixed value of [3 pi , the damping rate scales with y/T and this makes the damping strong in the case of galaxy 



II With this formula it is dW/dt = ~TW 
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clusters (T ~ 10 7 ~ 8 K). For (3 p i >> 1 the TTD damping rate from thermal electrons is T e /uj r « y^3irx/20 exp(— 5x/3) sin 2 8, 
where x = {m e /m p )/ cos 2 8, which is sufficiently small** to make the linear-theory approach adopted here still reasonable. 
Eq.(29) is a general expression of the damping rate due to TTD resonance with thermal particles from which well known 
formulae can be readily re-obtained. For instance in the case of low (3 p i it is V p h — > va and (l-B^p/ie-TrW 7 ) — > 1/2 and one gets 
the usual TTD-damping rate of fast modes with thermal electrons (e.g., Akhiezer et al. 1975; Achterberg 1981; Miller 1991): 

r m [Wm e vte Sm 2 8 

1 e/p 



[W m e vte sin 8 f v A \ . . 

p (k, 8) -> J | — exp <^ - — -= — }v A k (30) 

p \2m p VA \cos0\ I 2v 2 p cos 2 0j 



where we define vte = \J UbT /m e - A formula equal to Eq. 30 is also given for f3 p i « 1 in Ginzburg (1961) and Shafranov 
(1967) without adopting the simplified quasi-linear approach. These authors also report a non-quasi-linear formula for the 
damping rate of thermal electrons and protons under the particular condition of c s << V p h << v te , in which case the plasma 
dielectric tensor can be largely simplified by expanding the Z-function (Appendix B, Eqs. B17-B18) of electrons and protons 
for large (protons) and small (electrons) arguments. In this case the normalization of the formula for the damping rate of 
protons is 5 times larger than that in Eq. 30, while the formula for the damping of electrons is equal to Eq. 30 (this asymmetry 
in the electron-proton contribution comes from the expansions of the Z-function in the two opposite regimes for electrons 
and protons). Still since it is c s << V p h the damping from protons is negligible and Eq. 30 is equivalent to the result reported 
by these authors. 

The damping rate due to ultra-relativistic electrons and protons is given by (Appendix C, and BS73): 



r (kd\- |gfc|2 ( k± Yf k Wi Vph k \ K±/p V » h /,(i ( v * k V\ 2 n^riJ 9 ^ 



(31) 

e/p 



while the damping rate due to generic non-ultra relativistic and non-thermal particles is given in Appendix C (Eq.C9). 

In Fig. 4 we report the damping rate from both thermal and relativistic particles under conditions typical of massive 
and hot galaxy clusters. The most important damping for a mode propagating at small angles (8 < 1) is that with thermal 
protons, on the other hand, a mode propagating at larger angles is damped by thermal electrons. We find that under viable 
physical conditions the damping due to relativistic particles is formally relevant only in a narrow range of the values of 8 
(close to 8 ~ 7r/2), and that it accounts for only a few percent of the total damping rate. 

For a given temperature of the plasma, T, the strength of the damping rate decreases with decreasing j3 p i as the phase 
velocity of the modes increases with respect to the thermal velocity and this makes the particle-mode resonance more difficult. 

We notice that the overall damping rate is anisotropic with a relatively narrow peak at k/ku ~ 30 (for high j3 p i) where 
the bulk of thermal electrons resonates with the modes. On the other hand, as discussed in § 3 , the ICM turbulence is 
super-Alfvenic and thus the turbulent modes can easily bend the magnetic field lines. The time scale of the bending of lines 
from hydro-motions on a scale I is expected to be a fraction of I /vi , where vi is the rms velocity of the turbulent eddies 
at the scale I. The bending of the lines by hydro-motions on the shortest scales is thus the most efficient so that we can 
grossly estimate this bending time-scale, Tbb, as w a fraction of Ia/va; eddies on scales below I a cannot significantly bend 
the field lines. TT This value of Tbb should be compared with that of the damping time at collisionless scales which is grossly 
(from Eqs. 29 and 24) r(fc) _1 (/ m / p ) w \J m p /m e l m f v /c s - The relevant time-scale for isotropization of the pitch angle 8 due 
to line-bending is thus faster than the damping process, i.e. Tbb < l/r(i m /p), in the case " : 



>> 



mfp 



The condition in Eq.(32) is always satisfied in the ICM, at least under the hypothesis of this paper, and thus we shall use an 
effective damping rate for the bulk of the spectrum of magnetosonic modes which comes from the contribution from different 
8s and is defined by : 

/■f/2 

(T e/P (k))= / r e/p {k,8)sin8d8 (33) 
Jo 

This is reported in Fig. 5a as a function of c 2 /v A (— P P i/2) (for a given temperature of the ICM, see caption). Damping of 
magnetosonic modes is found to be always dominated by thermal electrons because they are faster than the phase velocity of 
these modes. The contribution from thermal protons drops for f3 p i < 20 since for smaller beta the phase velocity of the modes 
(Eq.4) increases with respect to the proton velocity and it is even more difficult for protons to satisfy the resonant condition. 



** Y e /u) r has a maximum value tti 0.2 

tt The wandering of the magnetic field at scales I < I a is discussed in Yan k, Lazarian (2004). 

it Here we assume that turbulent eddies reach scales < Ia (§ 5.1.3, Fig. 6a), in case l cu t > Ia the bending time-scale gets grossly of the 
order of a fraction of the damping time-scale at the cut-off scale, which would still be sufficient to have some isotropization 
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Figure 3. The ratio between particle energy and magnetic energy in the mode is reported as a function of 9; kgT=8.6 keV is assumed. 
Calculations are reported for: c 2 /v 2 ^=0.1, 1, 3 and 100 (from the bottom to the top of the diagram). The behavior at 9 ~ tt/2 is 
highlighted in the right panel for c 2 s /v 2 A = 100. 
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Figure 4. The damping rates of magnetosonic modes due to TTD resonance with thermal electrons (upper curves in the right end of 
the panel) and with thermal protons (upper curves in the left end of the panel), and the total damping rate (thick curves) are reported 
as a function of 6. The behavior at about 9 = w/2 is highlighted in the right panel. Calculations arc reported for : c^/i>^=100 (solid 
lines), and 1 (dashed lines), and taking k = 1 kpc -1 and ksT =8.6 keV. The damping rate with relativistic protons is also reported in 
both panels (dotted lines): in this case we assume an energy distribution in the form f(p) oc p~ 4 ' 2 and an energy density ~ 5% of the 
thermal one. 
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Figure 5. Panel a) The average damping rates of magnetosonic modes (Eq. 33), due to TTD resonance with thermal electrons (solid 
line) and thermal protons (dashed lines), and total damping rate (thick solid line) are reported as a function of c 2 /v 2 A {= (3 p i/2); k = 1 
kpc~ 1 is taken. Panel b) The turbulent cut-off scale (Eq. 45) is reported as a function of c 2 /v 2 A (= /3 p ;/2). Calculations are obtained 
assuming L o =300 kpc, and (Vl/c s ) 2 = 0.15 (upper curve) and (Vl/c s ) 2 = 0.3 (lower curve). In both panels we assume fcgT = 8.6 keV. 

Finally, let us comment that it is (T)/uj r << 1 and this further motivate the practical use of the quasi-linear theory in 
this paper. 

5 STOCHASTIC PARTICLE ACCELERATION IN GALAXY CLUSTERS 

In this Section we discuss the particle acceleration process in the ICM via resonant and non-resonant mechanisms with 
compressible modes. 



5.1 Resonant Transit Time Damping Acceleration 

5.1.1 Introduction 

Compressible (and incompressible) low-frequency MHD waves can strongly affect particle motion through the action of the 
mode-electric field via gyroresonant interaction (e.g., Melrose 1968), the condition for which is : 

u) — fc||V|| — n— = (34) 

where n = ±1, ±2, .. gives the first (fundamental), second, .. harmonics of the resonance, while v\\ = fj,v and fcy = nk are 
the parallel (projected along B ) speed of the particles and the wave-number, respectively. In general gyroresonance is a 
process important only for modes at very small scales, I « I a- However, as anticipated in § 3.5, at these scales fast modes 
are probably absent in the ICM due to strong resonant dampings (§ 5.1.3, Figs. 4, 5, 6) and because they do not couple with 
the Alfvenic cascade (§ 3.3.2, 3.5). 

Interestingly enough, the compressible component of the magnetic field of compressible modes (i.e. the component along 
B in the case of oblique propagation) can interact with particles through the n = resonance. This interaction is called 
transit-time damping (e.g., Fisk 1976; Eilek 1979; Miller, Larosa & Moore 1996; Schlickeiser & Miller 1998). An important 
aspect of this interaction is the need of isotropization of particle momenta during acceleration (e.g., Schlickeiser & Miller 1998). 
This is because the n = resonance changes only the component of the particle momentum parallel to the seed magnetic 
field. This would cause an increasing degree of anisotropy of the particle distribution and thus the deriving acceleration would 
become less and less efficient with time. Under our working picture, particle-pitch angle scattering in the ICM can be provided 
by several processes discussed in the literature. Those include electron firehose instability which is indeed driven by pressure- 
anisotropies in high beta plasma (Pilipp & Volk 1971; Paesold & Benz 1999), and gyro- resonance by Alfven (and slow) modes 
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at small scales, provided that these modes are not too much anisotropic (cf. Yan & Lazarian 2004). The latter condition means 
that the Alfvenic modes are considered for scales not much less than I a, provided that the turbulence injection is isotropic. In 
addition, gyro-resonance was discussed for the electrostatic lower hybrid modes generated by anomalous Doppler resonance 
instability due to pitch angle anisotropies (e.g., Liu & Mok 1977; Moghaddam-Taaheri et al. 1985) and by gyroresonant 
interaction with whistlers (e.g., Steinacker & Miller 1992). The latter process, however, is somewhat more problematic than 
the Alfvenic mode scattering, as whistler turbulence is even more anisotropic than the Alfvenic one (Cho & Lazarian 2004). 
Finally, instabilities within cosmic ray fluid look as a safer bet for isotropizing cosmic rays. For instance, Lazarian & Beresnyak 
(2006) proposed isotropization of cosmic rays due to gyroresonance instability that arises as the distribution of cosmic rays 
gets anisotropic in phase space. This instability that is customary discussed for plasma rather than for cosmic rays (see Gary 
et al. 1994, Kulsrud 2004) would guarantee that in the environments of galaxy clusters the TTD will not be quenched. 



5.1.2 Diffusion Coefficient 

The momentum-diffusion coefficient, D pp , of particles can be calculated by deriving the first-order corrections due to small 
amplitude plasma turbulence to the orbits of particles in a uniform magnetic field, and ensemble averaging over the statistical 
properties of the turbulence (e.g., Jokipii 1966). The resulting analytic expressions for the pitch-angle and momentum- 
diffusion coefficients due to TTD resonance with fast modes in a low beta plasma can be found in Schlickeiser & Miller (1998). 

An additional and self-consistent way to derive the momentum-diffusion coefficient from the quasi-linear theory is to 
use an argument of detailed balancing. The diffusion coefficient of a a-species is indeed related to the damping rate of the 
modes themselves with the same particles, and one has (e.g., Eilek 1979; Achterberg 1981) : 

J d 3 P E a (^f^J = J dkr a (k,6)W(k) (35) 

where E a is the energy of a particle of species a, and W(k) is the total energy of the modes in the elemental range dk. This 
is given by : 

where (W/We)^ is the ratio between the total and the electric energy in a single mode propagating at k (§ 4.2), and We(/c) 
is the electric-field energy of the modes in the elemental range dk. In Eq.(36) we have assumed an isotropic spectrum of 
the electric field fluctuations which is an appropriate assumption for super-Alfvenic turbulence and fast modes (e.g., Cho & 
Lazarian 2003). 

If isotropy of the particle momenta is maintained, the time evolution of the particle distribution function is related to the 
diffusion coefficient by : 

OUP) _ 1 ( p * D JjM) (37) 



dt p 2 dp y dp 

and thus Eq.(35) reads 



/ 



" (v'^-W 1 ) -i/*/ *»«-)r«( M)W «, t ) ( £) (m) (38) 



Here we are interested in deriving the diffusion coefficient in the case of relativistic species in the ICM (§ 5.1.3, 5.5). The 
damping with these particles (Eq. 31) can be expressed in the form : 

F a (k, 6) = -T a (6) k J p 4 dp^0^- (39) 
and from partial integration of Eq.(38) and from Eqs.(39) and (31) taking Ws(fc) = (c/V p h) 2 WE(k), one gets : 



r/2 



This represents a self-consistent average (in terms of particle pitch-angle) momentum-diffusion coefficient of isotropic particles 
with momentum p which couple with fast magnetosonic modes via TTD resonance. Eq.(40) in its low beta plasma limit 
(essentially Vph — > va and V p h « c) is consistent with the expression (Eq.29) given in Schlickeiser & Miller (1998) in its 
z — k±v±/fl << 1 limit and averaged over the particle pitch-angle §§. 

It is sufficient to integrate (average) Eq.(29) in Schlickeiser & Miller (1998) over the particle pitch-angle using the properties of the 
delta-function, to solve this integration and to expand the Bessel function in Eq.(29) for small arguments. 



(0)^(, V ph /c\ ( .V ph /c 
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5.1.3 Acceleration efficiency in the ICM 

As summarized in § 3.5 we focus on a picture in which compressible turbulence is injected at large scales by the action 
of cluster mergers and accretion of matter. Provided that large scale turbulence in the ICM is not significantly affected by 
the ion-viscosity (§ 3.4), an inertial range is established due to the combination of turbulence injection and cascading. For 
isotropic turbulence the diffusion equation in the k-space is given by : 

^ ' i 

where D kk is the diffusion coefficient in the k-space, r t (k,t) are the different damping terms (§ 4.3), and I(k,t) accounts 
for the turbulence injection term. The wave-wave diffusion coefficient of magnetosonic modes (Kraichnan treatment; see also 
Zhou & Matthaeus 1990; Miller, La Rosa, & Moore 1996 for low beta plasma) is given by 

We assume a constant (in time) injection spectrum of the modes in the simple form I(k) — I 8(k — k a ) so that the stationary 
spectrum of turbulence at the scales not significantly affected by dampings (r, ~ 0) can be readily obtained from Eq.(41) : 

W(k)= (p oP (V ph )y k~i (43) 
and the cascading time at a the scale I = 2ir/k, is given by : 
2 ^7 <V,k>A * 



Provided that the dissipation of compressible turbulence in the ICM is collisionless (§ 3.4), the turbulence cascading gets 
suppressed at a scale at which the resonant damping time-scale, P -1 , approach the cascading time. This scale is given by 
Eq.(44): 



81 Iq ( {T{k)) X 
Up{V ph ) ^ fc J 



(45) 



where (P) is the average collisionless TTD damping term given by Eqs.(29), (31) and (33). The value of the cut-off scale is 
reported in Fig. 5b as a function of the beta of the plasma for physical conditions in the ICM (see caption): we find that if 
turbulence is energetic enough (actually for the values used in the § 5.3-5.5) compressible modes are dissipated at « sub-kpc 
scales. The cut-off scale slightly increases in the case of small (3 p i as the cascading of magnetosonic modes becomes less efficient 
(the cascading time-scale goes as T k k oc (Vph) and, fixed c s , increases for small f3 p i). Actually the cascading of compressible 
motions is likely to reach MHD scales before being dissipated, l cut < I a (Fig. 6a), and in this case it is also worth to mention 
that an Alfvenic turbulence can be activated by the cascading of these compressible motions. 

Eq.(41) is appropriate to describe the time evolution of the total spectrum of isotropic turbulent modes. On the other 
hand, formally in the collisionless regime the ratio between the energy of the fields (E and B) and that associated with 
particles changes with the mode-propagation angle (Figs. 1-2, Appendix B). However the induced anisotropy is within a 
factor of 2-3 for a stationary B , and it should be efficiently smoothed out by the effect of the bending of the field lines 
(§ 4.3). Thus we shall adopt isotropy as a viable approximation, and define the energy associated with the magnetic field 
fluctuations as: 

WB ^~^^^) mk ' t] (46) 

where for consistency |Bfc| /W is taken from Eq.(24) and () indicates the average over the propagation angle of the modes. 
The TTD diffusion coefficient in the particle-momentum space is then obtained from Eqs.(40), (43), and (46) in the form : 

np\ MB k \\ 1 ( 2I {V ph ) \ l2 k (f) | r^^tJi V ph /c\ ( f Vph/c 2 y 

DvAp ^ - z-^^Mw^K^p— ) kcut[t) J d9Vph j^m n y 1 y 1 { ^e> ) (47) 

Eq.(47) allows a prompt estimate of the acceleration efficiency via TTD resonance, once the injection rate per unit mass of 
the compressible turbulence (I /p) and the injection scale, k (or L ), are fixed : 



^ Here (Vph) is essentially a representative, averaged (with respect to 9) phase velocity. 
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where C « 5 — 6 is a numerical factor which can be readily obtained by taking I a /p w VI/tll and Eq.(44). The resulting 
systematic acceleration rate, T aC c, is given by : 

Tacc= P 3 {-?-(p 2 D pp (p))}- 1 (49) 

The systematic acceleration time from TTD resonance does not depend on particle momentum (see also Fig.7) and is reported 
in Fig. (6b) as a function of c 2 /v 2 i (= 2/3 p i): for a given temperature (and f3 p i > 1) the acceleration efficiency scales approxima- 
tively with y/T and is found to be almost independent from the value of (3 p i . The important point here is that the strength of 
the TTD-acceleration efficiency, powered by compressible turbulence with large-scale rms velocity V L /c a « 0.3, is found to 
give a systematic acceleration time of the order of ~ 10 8 yrs which is sufficient to accelerate electrons up to energies of several 
GeV, and this may produce diffuse synchortron radio emission in /iG-magnetized media (§ 6). 



5.2 Nonresonant acceleration 

5. 2. 1 Introduction 

Resonant TTD acceleration is not the only process by which compressible turbulence may accelerate cosmic rays in the ICM. 
For instance, fast particles can be accelerated also by large scale compressible motions (e.g., Ptuskin 1988; Chandran 2003; 
Chandran & Maron 2004; Cho & Lazarian 2006) . Compression changes the particle momentum according to : 

th-JpV.V, (50) 

If the medium is neither expanding nor contracting it is (V ■ V;) = and thus particles will not experience regular changes in 
energy. On the other hand if V; is a turbulent field a statistical acceleration effect (analogous to a classical second order Fermi 
process) may exist. This is essentially because particles would statistically experience more compression than expansion. 



5.2.2 Diffusion Coefficient 

Limiting to the case V 2 « c 2 s and provided that the turbulent velocity of the medium has correlation scales much longer than 
the effective particle mean free path, the diffusion coefficient in the particle momentum space, D pp , and the total (turbulent 
advection and diffusion) spatial diffusion coefficient, D*, can be obtained by standard procedure in plasma physics in the 
quasi-linear approximation. These are (Ptuskin 1988): 



dyy 2 V(y) 



and 



^=D[l + tfjm,) (52) 



3 J k c 2 + y*D» / 

where D is the spatial particle-diffusion coefficient (without considering the effects induced by the nonresonant compressible 
coupling itself, Eq. 51), and V(y) is defined as : 



/' 



V(y)dy = Vt (53) 

In this regime slow and fast diffusion limit exist. In the slow limit the rate of particle diffusion out of compressible eddies is 
slower than the wave period, t w ~ l/c s , i.e. Tdiff ~ l 2 /D >> t w and c 2 >> k 2 D 2 . Here the process is mainly contributed by 
the action of the smaller eddies in the spectrum of the modes and it becomes faster as this minimum scale gets smaller (e.g., 
Cho & Lazarian 2006). From Eq.(51) we find : 

For small minimum-turbulent scales this process formally becomes extremely efficient, however the minimum scale of the bulk 
of compressible turbulent eddies in the ICM cannot be very small as these modes are strongly damped (§ 5.1.3, Fig. 5). 
In the opposite case, in the fast diffusion limit, particles leave the eddies before they turnover, i.e. Tdiff << t w and c 2 « k 2 D 2 . 
Here the process is mainly contributed by the action of the largest eddies which contain the bulk of the turbulent energy, and 
from Eqs.(51) & (53) we find : 

D PP ~ \p 2 ^ (55) 
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Figure 6. Panel a): The ratio lA/lcut is reported as a function of B. Calculations are reported for fcgT =6 keV (dashed lines), 9 keV 
(dotted lines), and 12 keV (solid lines); (Vl/c 3 ) 2 = 0.3, n t h = 10 _3 cm — 3 , and L = 300 kpc were assumed in the calculations. 
Panel b): The acceleration time (Eq. 49) is reported as a function of c^/v 2 ^ in the case of nonresonant compressive acceleration 
(Eqs. 51) and resonant TTD acceleration (Eq. 47). Calculations are reported for kgT =7 keV (dashed lines) and kgT =11 keV (solid 
lines); (V L /c s ) 2 = 0.3, n th = 10~ 3 cm~ 3 , and L a = 300 kpc were assumed in the calculations. The acceleration time from the combined 
effect of the two mechanisms is also shown (tick lines). 

The box marks the relevant range of the values of c 2 /v 2 A {= (3 p i/2) in the hot ICM, and the acceleration time necessary to boost relativistic 
electrons at several GeV (this accounts for both synchrotron and inverse Compton losses with redshift ~0-0.3). 



An important point discussed in § 3.2 is that particle-spatial diffusion itself is likely to be affected by the turbulent bending 
of the magnetic field lines which gets the effective ion mean free path ~ I a- Compared to the Coulomb or gyroresonance 
scattering the diffusion with the characteristic scale Ia does not involve any changes of the particle energy via scattering. 
Therefore the particle may diffuse slowly, but the only change in energy will be due to large scale compressions (cf. Cho & 
Lazarian 2006). We thus shall adopt a very simplified form of the spatial diffusion coefficient in Eqs. (51-55) : 

D ~ |/3max |/ cut , min{/A , L/p}} (56) 

The combination between Eqs. (51) and a turbulent -driven spatial diffusion coefficient (e.g., Eq. 56) provides an important 
refinement of the evaluation of the cosmic-ray acceleration via compressible long wave turbulence, and may have important 
consequences in the case of the particle acceleration in the ICM (§ 5.3). 

Finally, we want to remind that Eq.(51) is obtained by neglecting the effect of possible additional scattering processes due to 
resonant particle-wave interactions. The presence of instabilities in cosmic rays may create an additional slab-type Alfvenic 
component that would produce additional gyroresonance acceleration and reduce the effective mean free path (Lazarian & 
Beresnyak 2006). Conservatively we do not discuss this possibility in the present paper. 



5.2.3 Acceleration efficiency in the ICM 

In this Section we calculate the efficiency of the particle acceleration from large-scale nonresonant compression in the ICM. 
Taking a Kraichnan scaling for the super-Alfvenic compressible turbulence, V(fc) « V 2 k~^l 2 / L^J 2 , from Eq.(51) we have : 



D s Dp iPU (57) 

where the spatial-diffusion coefficient is given by Eq.(56). 

The resulting systematic acceleration time is independent of particle momentum (at least in the ultra-relativistic case, see 
also Fig. 7) and is reported in Fig. 6b. For a given temperature of the plasma, T, in the case of small /3 p i the nonresonant 
compression is formally very inefficient because for large values of the magnetic field the particle spatial-diffusion coefficient 
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is large (essentially D « 1/3 (3cl m f p , l m f P from Eq. 1). On the other hand, in the case of large [3 p i the acceleration efficiency 
increases because turbulence bends the magnetic field lines at scales smaller than l m f p and the effective particle mean free 
path is w I a (which scales as P~^ 2 )', saturation for large (3 p i is reached when l cut > I a (Fig- 6). 

The reference value of (3 p i in the ICM is in the range f3 pl ~ 200 - 1000 (i.e. B ~ 0.5 - 3/uG, with n th ~ 10~ 4 - 10~ 3 cm~ 3 and 
fcsT ~ 7 — 10 keV), and formally under these conditions we find that the acceleration efficiency from nonresonant compression 
driven by relatively energetic turbulence (caption) is similar to that due to the TTD-resonance. 

As already pointed out in § 5.2.2, in the derivation of Eq.(51) (or Eq.57) it was assumed that the effective particle mean free 
path is much smaller than the scale of the turbulent eddies. This condition is formally violated in the case of small j3 p i in 
Fig. 6b where the smaller turbulent eddies are < l m fp (mean free path l m f P ~ 10 — 50 kpc). On the other hand, this does 
not happen for larger f3 p i, since in this case the particle effective mean free path, « I a ~ hut, is actually comparable to (or 
smaller than) the smallest turbulent eddies. 

5.3 Overall effect of compressible turbulence 

As discussed in § 5.1.1 the TTD resonance is expected to be an efficient mechanism in the ICM, provided that particle 
isotropy is preserved. Yet the TTD alone might not be efficient enough in maintaining such isotropy because both D pp (fi,p) 
and D MM (/x,p) are strongly maximized for particles moving at small angles with the direction of the seed magnetic field. 
However additional resonant processes acting on small scales might easily maintain particle isotropy. If these mechanisms 
are really at work in the ICM they should also affect the spatial diffusion, D, of the particles and thus the efficiency of the 
nonresonant compression mechanism. Formally with decreasing D the nonresonant coupling with eddies in the fast diffusion 
limit becomes more efficient, and, at the same time, a larger range of scales of the eddies couples with particles in the slow 
diffusion regime which is very efficient; actually this is what happens with increasing the beta of the plasma in Fig. 6b. 
However, if the spatial diffusion is strongly suppressed, namely when D < c s l cut in Eqs.(51) and (57), one gets into the slow 
diffusion limit at any turbulent scale, and a decrease of D yields a corresponding decrease in the efficiency of the nonresonant 
compression (Eq. 54). Thus future studies using self-consistent spatial diffusion coefficients will be of great importance. 

The turbulent bending of the field lines which happens in the super-Alfvenic case cannot change the pitch angle of particles 
which would preserve the adiabatic invariant, however in the high beta ICM turbulent bending is associated with turbulent 
compressions which indeed power the nonresonant acceleration mechanism and might provide a source of particle-pitch angle 
isotropization. The spatial diffusion coefficient is related to that in the pitch angle as (order of magnitude) D w c 2 /D MM , 
and the resulting time-scale of the pitch angle scattering, w D^, is indeed much shorter than the acceleration time of fast 
particles (which is w 10 7 — 10 8 yrs). 

This is important since it implies that the action of large scale compressible turbulence in the ICM is twofold. On one hand 
particles diffusing through the compressible turbulent eddies experience substantial nonresonant stochastic acceleration. On 
the other hand, even without requiring additional processes at small scales, this might contribute to help in maintaining 
particle-momentum isotropization, so that the compressive component of the turbulent magnetic field (that along B ) may 
also couple efficiently with particles via TTD resonance without greatly change the particle spatial diffusion. 
These two mechanisms, TTD resonance and nonresonant compression, are driven by the same turbulent modes and involve 
independent particle-mode couplings and thus, as a first approximation, the acceleration process may be thought as the 
combination of the two effects; the deriving particle acceleration time is also reported in Fig. (6b) . 



6 COMPRESSIVE TURBULENCE AND PARTICLE RE-ACCELERATION MODEL IN GALAXY 
CLUSTERS 

As already anticipated in the Introduction direct evidence for relativistic electrons diffused on Mpc scales in the ICM comes 
from radio halos and relics (e.g., Feretti 2005), while the hard X-ray tails detected in a few clusters may result from inverse 
Compton scattering of the Cosmic Microwave Background photons by the same electrons (e.g., Fusco-Femiano et al. 2004; 
Rephaeli, Gruber & Arieli 2006). 

The particle re-acceleration model is a promising possibility to explain the properties of the giant radio halos and possibly 
also the strength of the hard X-ray tails. This scenario assumes that turbulence is injected in a substantial fraction, Mpc 3 , of 
the cluster volume during cluster-cluster mergers, and that relativistic electrons already present in the ICM and accumulated 
at 7 w 100 are re-accelerated for a typical time-scale of <Gyr (e.g., Brunetti et al. 2001,04; Petrosian 2001; Fujita et al. 2003). 
Alternatively these seeds electrons to be re-accelerated could be secondary products of hadronic interactions (Brunetti & Blasi 
2005). 

In this Section, after a brief review of the injection processes of cosmic rays in galaxy clusters and of the most relevant 
channels of energy losses (§ 6.1), we provide calculations in the context of the particle re-acceleration model which include 
the effect of TTD-resonance and nonresonant acceleration due to compressible turbulent modes injected at large scales. 
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6.1 Cosmic Ray injection in the ICM 

There is a general consensus on the fact that several mechanisms of injection of cosmic rays may be at work in the ICM, and 
that once injected the bulk of these cosmic rays does not escape the cluster (e.g., Berezinsky, Blasi & Ptuskin 1997; Ensslin 
et al. 1998; Voelk & Atoyan 1999). 

Collisionless shocks are generally recognized as efficient particle accelerators through the so-called diffusive shock acceler- 
ation (DSA) process (Drury, 1983; Blandford & Eichler 1987). This mechanism has been invoked several times as an efficient 
acceleration process in clusters of galaxies (Takizawa & Naito 2000; Blasi 2001; Miniati et al. 2001; Fujita & Sarazin 2001; 
Ryu et al. 2003). Present simulations confirm the analytical claim that shocks with Mach number larger than 2-3 are rare 
(Gabici & Blasi 2003), and claim that the energy content in the form of cosmic rays in massive clusters may be of the order 
of a few percent of the thermal energy (Pfrommer et al. 2006; Jubelgas et al. 2006). The bulk of the energy of these cosmic 
rays is injected in the cluster outskirts by shocks with a Mach number of the order of ~ 3, the real efficiency of these shocks is 
however uncertain and it is generally computed in the test particle limit and according to the so-called thermal leakage model 
(e.g., Kang & Jones, 1995). 

A contribution to the injection of cosmic rays in clusters of galaxies may come from Active Galactic Nuclei which indeed 
might fill the ICM with relativistic particles and magnetic fields, extracted from the accretion power of their central black hole 
(Ensslin et al., 1997). Similarly to Active Galactic Nuclei, powerful Galactic Winds may also inject relativistic particles and 
magnetic fields in the ICM (Volk & Atoyan 1999). Although the present day level of starburst activity is low, it is expected 
that these winds were more powerful during starburst activity in early galaxies, as also suggested by the iron abundances in 
galaxy clusters (Volk et al. 1996). 



6.2 Energy Losses 

6.2.1 Electrons and Positrons 

In the conditions typical of the ICM, ultra-relativistic electrons rapidly cool down through inverse Compton and synchrotron 
emission, and accumulate at Lorentz factors 7 ~ 100 — 500 where they may survive for a few billion years before cooling further 
down in energy through Coulomb scattering and eventually thermalize. Energy losses and relevant time-scales of relativistic 
electrons in the ICM are discussed in several papers (e.g., Sarazin 1999; Petrosian 2001; Brunetti et al. 2004; Pfrommer & 
Ensslin 2004). In Fig. 7a we report the particle life-time as a function of the Lorentz factor: the life-time has a peak at 
7 ^ 10 2 — 10 3 where the cooling of electrons is slower and where particles may accumulate providing a seed populations to be 
re-accelerated in the context of the re-acceleration model. More specifically, Fig. 7a is obtained for typical physical conditions 
in cluster cores and in the cluster outskirts: in the external regions of clusters electrons survive since Coulomb losses are less 
severe and in principle these particles can be accumulated for cosmological time-scales at energies 7 ~ 100 — 1000. On the 
other hand, in cluster cores the higher thermal density limits the maximum life-time of electrons at less than 1 Gyr. 



6.2.2 Protons 

Once injected the relativistic cosmic-ray protons do not suffer catastrophic radiative-energy losses. The only relevant channel 
of energy losses for these particles in the ICM is given by hadronic collisions which however get a typical particle life-time 
which is larger than a Hubble time for ~ GeV particles. This, together with the long time necessary to the bulk of these 
cosmic rays to diffuse out of clusters, makes clusters themselves reservoir in which cosmic ray protons are confined and may 
accumulate over cosmological epochs (e.g., V61k et al. 1996; Berezinsky, Blasi & Ptuskin 1997). 

On the other hand, mildly and sub- relativistic protons may be significantly affected by Coulomb energy losses, which in turn 
change the particle spectrum with respect to the injection spectrum. The rate of Coulomb losses is (e.g., Schlickeiser 2002) : 



r°'*«"(-» 2 )-i( i+ s)« [ '(- < i )2 > 



(58) 



where f3 p and f3 e ~ 0.18(T/10 s i(') 1 ^ 2 are the velocity in units of the light speed of thermal electrons in the ICM and of the 
cosmic ray protons, respectively. 

As in the case of leptons, the details of the mechanisms of energy losses of cosmic ray hadrons in the ICM can be found in 
several papers (e.g., Blasi & Colafrancesco 1999; Pfrommer & Ensslin 2004; Brunetti & Blasi 2005). In Fig. 7b we report the 
particle life-time as a function of the particle momentum. Fig. 7b is obtained for typical physical conditions in cluster cores 
and in the cluster outskirts: it is clear that even in the cluster cores where losses are much severe, the bulk of relativistic 
protons has a life-time of the order of an Hubble time. Only protons with kinetic energy larger than about 200 GeV and 
smaller than about 30 MeV in the cluster cores have life-times smaller than a couple of Gyrs, while just out of the core regions 
the life-time of these particles grows (time oc n^ 1 ) and all these particles are expected to survive for cosmological time-scales. 
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Figure 7. Left panel: The life— time of relativistic electrons in the ICM at z = 0.2 as a function of the Lorentz factor. Thick (blue) lines 
are for cluster cores and thin (red) lines are for cluster periphery. We report the total life-time (solid lines) and the life— times due to 
single processes: Coulomb losses (dashed lines), synchrotron and IC losses (dotted lines), and bremsstrahlung losses (long dashed lines). 
Right panel: The life— time of cosmic— ray protons in the ICM at z = 0.2 as a function of the particle momentum. Thick (blue) lines are 
for cluster cores and thin (red) lines are for cluster periphery. We report the total life— time (solid lines) and the life-times due to single 
processes: Coulomb losses (dotted lines) and pp-collisions (dashed lines). 

In both panels calculations in the cores (thick— blue) are obtained for B = 3(j,G and n th = 2 X 10~ 3 cm~ 3 , and in the periphery (thin— red) 
for B = 0.5/iG and n t h = 10 — 4 cm — 3 . For comparison, the dash-dotted lines in both panels give the acceleration time— scale which is 
used in Fig. 8; note that the increase of this time-scale in the case of sub— relativistic protons is due to the decrease of the efficiency of 
the non— resonant compression at sub— relativistic energies. 



6.3 Numerical Calculations 

In this Section we calculate the time-evolution of the spectrum of the relativistic particles stochastically re-accelerated by 
turbulent modes in the ICM. 



6. 3. 1 Formalism 



As discussed in § 5 we shall assume isotropy of the particle momenta and of the modes, and in this case the time evolution of 
the spectrum of the turbulent modes and of the particles can be formally derived by a set of coupled kinetic equations. The 
time evolution of the spectrum of the leptonic component is given by : 
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(59) 



where N- and N+ stands for electrons and positrons (JV = 47rp 2 /, / is used in § 4), respectively, and where the terms 
dp/dt ra d and dp/dti account for radiative (synchrotron and IC) and Coulomb losses, while D pp and D pp are the resonant 
(TTD, Eq.47) and the non-resonant (from turbulent -compression, Eq. 57) particle momentum-diffusion coefficients. In Eq. 59 
we also formally include an injection term, Q f ,±(N p (p,t);p,t), which depends on the spectrum of cosmic ray protons, and 
is necessary in the case that the re-accelerated electrons and positrons are injected by hadronic collisions in the ICM (see 
Brunetti & Blasi 2005). 

The time evolution of the spectrum of cosmic ray protons is given by : 
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(60) 



where the term dp/dti accounts for Coulomb losses (Eq. 58), while the particle depletion due to proton-proton collisions can be 
neglected. Because the population of cosmic ray protons in the ICM essentially comes from the accumulation of these particles 
during cosmological time-scales, we do not consider the source term in Eq.(60) which would account for the contribution from 
freshly injected protons. 
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Figure 8. Left Panel: Time-evolution of the spectrum of relativistic electrons as a function of the Lorentz factor. 
Right Panel: Time-evolution of the spectrum of cosmic ray protons as a function of the particle momentum. 

In both panels calculations are reported for: t = 0, 4 X 10 15 , 8 X 10 15 , 10 16 , 1.2 X 10 16 sec from the start of the re— acceleration phase. 
Calculations are performed assuming (V L /c s ) 2 = 0.18, L = 300 kpc, n th = 10~ 3 , k B T=9 keV, B = 1/j.G, and redshift z=0.1 (for IC 
losses) . 

The evolution of the spectrum of the compressible turbulent modes is given by Eq.(41) described in § 5.1.3, where all the 
dampings are formally derived in combination with Eqs. (59-60). The effect of the non-resonant damping on the spectrum of 
the modes can be neglected since turbulent-compression acts efficiently only on relativistic particles in the ICM (§ 5.2), and 
this gets a net damping rate which is much smaller than that via TTD-resonance with the thermal ICM. 



6.3.2 Assumptions 

In this paper we adopt the particle re-acceleration model assuming that a seed population of relativistic electrons and protons 
in the ICM is re-accelerated by turbulence injected at large scales during a merger event. For simplicity we do not study 
the more complex issue of the re-acceleration of secondary electrons and positrons injected by proton collisions in the ICM 
(Brunetti & Blasi 2005). 

The new of this paper is that compressible turbulence is used as the driving of particle re-acceleration, and accordingly the 
detailed diffusion coefficients obtained in § 5 and the scenario and properties of turbulence discussed in § 3-4 are used in the 
calculations. For seek of clarity the main assumptions and physical parameters used in the calculations are listed below : 

i) We consider physical parameters appropriate for massive galaxy clusters: T w 10 8 K, n t h ~ 10~ 3 cm~ 3 , B w 0.5 — 3 fiG. 

ii) Turbulence is assumed to be sub-sonic, with V£ « c 2 s , and injected at large scales Lo ~ 300 — 500 kpc for a typical 
cluster-cluster crossing time. 

ill) The initial spectrum of electrons, N e (p,t = 0), is derived by assuming that electrons are injected in the ICM in a 
single event and then evolve passively for w 1-3 Gyr before being re-accelerated (see Brunetti et al. 2004). 

iv) The initial spectrum of protons, N p (p,t — 0), is derived by assuming that protons are continuously injected in the 
ICM for a long period, ^3 — 5 Gyr, with a constant injection spectral rate Q oc p~ 2 ' 2 before being re-accelerated (see Brunetti 
et al. 2004). 

v) Damping terms from relativistic species in Eq.(41) are neglected in the calculations as they become important only if 
the relativistic component gets a relevant fraction of the thermal energy of the ICM (§ 4.3). 

vi) The damping term due to thermal particles is taken stationary because, under the assumption (ii), the thermal 
properties of the ICM are not significantly modified with time. 

Under conditions ii), v) and vi) the spectrum of the modes is also stationary and this is given by Eq.(43). 
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6.3.3 Main Results 

Once large scale turbulence is injected in the ICM, magnetosonic modes take a relatively long time to cascade at collisionless 
scales : 



In the re-acceleration scenario this is an unavoidable temporal gap, of a fraction of a Gyr, between the injection of the first 
turbulent eddies and the beginning of the particle re-acceleration process. When turbulence reaches collisionless scales the 
acceleration process starts and particles take a time, of the order of the re—acceleration time, to be significantly boosted in 
energy. A relevant example of the time evolution of the electron and proton spectrum during the re— acceleration period is 
reported in Fig. 8 assuming V% ~ 0.18 c 2 (see caption): the seed electrons initially accumulated at 7 ~ 10 2 — 10 3 are efficiently 
re-accelerated up to 7 « 10 4 — 10 5 . 

Radio (and hard X-rays) observations can be well explained in terms of a high energy tail of emitting relativistic electrons 
at energies of several GeV (e.g. Schlickeiser et al. 1987; Brunetti et al. 2001; Petrosian 2001). We calculate the evolution of 
the re-accelerated particles under the conditions given in § 6.3.2 and find that, quite independently from the initial electron 
and proton spectrum, an appreciable high energy tail of relativistic electrons at these energies is produced approximatively 
for (Vl/c s ) 2 > 0.15(1 + 2/O.I) 2 II II. In this case compressible turbulence injected at large scales in galaxy clusters may actually 
trigger efficient particle re-acceleration, and potentially explain the diffuse Mpc radio sources observed in massive galaxy 
clusters and the hard X-rays in excess to the thermal X-ray emission. 

A spectral break, at ^GHz frequencies, is observed in the synchrotron spectrum of a few radio halos and this is interpreted in 
favour of the re-acceleration scenario (e.g., Brunetti 2004; Feretti 2005). Under the conditions i) in Sect. 6.3.2, such a break 
requires the presence of a corresponding break in the spectrum of the emitting electrons at energies «5-10 GeV, and we find 
that this is reproduced by our re-acceleration model in the case of moderate turbulence, typically (Vl/c s ) 2 ~ 0.15 — 0.25, while 
in the case of more energetic turbulence electrons can be re-accelerated at larger energies and the corresponding synchrotron 
break is shifted at several GHz. 

Also protons are efficiently re-accelerated. Relativistic protons are not subject to radiative losses and since the re-acceleration 
efficiency scales with the energy of the particles (§ 5.1-5.2) the spectrum is simply shifted at higher energies and the slope of 
the injection spectrum is essentially preserved during the re-acceleration (Fig. 8). 

Under our assumptions (v), § 6.3.2), it is Y t h » Y rei and the spectrum of the turbulent fluctuations in terms of magnetic 
field, Wb , does not depend on the presence of cosmic rays, thus protons cannot significantly affect the acceleration process of 
relativistic electrons. This marks an important difference with Alfvenic re-acceleration, in which case the dominant damping 
of the modes comes from the resonance with relativistic protons and thus these protons affect the electron acceleration 
(Wave-proton Boiler, Brunetti et al. 2004). 

An additional point to stress here is that, because T t h >> Y re i, the fraction of the turbulent energy which goes into the 
cosmic rays via TTD-resonance is simply « Y re i/Y t h and is fixed by the fraction of the energy in the ICM which is in the 
form of cosmic rays. An additional contribution to the energy of the re-accelerated particles comes from the non-resonant 
compression. In our calculations (assuming that cosmic rays store a few percent of the thermal energy in the ICM) the total 
fraction of the turbulent energy which goes into non-thermal particles is of the order of w 2 — 5%. 

In case of long re-acceleration periods, actually > 3 — 4 times the re-acceleration time (Eq. 49), a non negligible fraction 
of the electron number is boosted toward the maximum energy. Here an equilibrium between acceleration and losses is reached 
and most of the energy flux from the damping of the turbulence with these particles is radiated away via synchrotron and IC 
by the same re-accelerated particles. Thus in principle for very long re-acceleration periods the total energy of the electron 
population should saturate and the spectrum is expected to slowly approach stationary conditions. On the other hand since 
cosmic ray protons are free from energy losses, the energy flux from the damping of the turbulence is totally stored in the form 
of particle energy and this gives an unbalance between electron and proton acceleration. In the re-acceleration scenario this 
unbalance is not expected to be large, indeed turbulence is injected during cluster mergers and the duration of a re-acceleration 
period is constrained by the cluster-cluster crossing time and by the turbulence cascading time, and these cannot significantly 
exceed about 1 Gyr (see Eq. 61). In addition, present studies of the number counts of giant radio halos in galaxy clusters 
limit the life-time of these sources at about « 1 Gyr (e.g., Hwang 2004) and this additionally constraints the duration of 
stochastic particle re-acceleration periods in galaxy clusters. Actually, given these limits, we find that assuming (Vl/c s ) 2 ~ 
0.15-0.25 (which is required to provide the necessary electron re-acceleration up to « 5 — 10 GeV) and a duration of the 
re-acceleration phase in the range 0.4-1 Gyr, the total energy of the cosmic ray protons and that of the relativistic electrons 
are both boosted by a factor 1.5-4, and the unbalance is not substantial. 




(61) 



II The term (1 + z/0.1) comes from IC losses (as for B /xG synchrotron losses are sub-dominant). 
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7 DISCUSSION 
7.1 Major Results 

The problem of proton and electron stochastic re-acceleration by compressible motions is a complex one. The efficiency of 
acceleration depends on the spectrum of compressible turbulent motions. The extend and the shape of this spectrum, in its 
turn depends on the processes of plasma damping. In the case of the hot ICM the corresponding issues have not yet been 
clarified sufficiently in the literature. 

As a result, we had to address those issues one by one. Namely, we started with the problem of describing turbulence in 
ICM in § 3. First of all, we provided arguments suggesting that turbulence is expected to be present in the medium in hot 
(and massive) galaxy clusters. This is also because the ICM is magnetized and this implies a partial suppression of the plasma 
viscosity. The suppression of the viscosity in a magnetized medium is a well known effect and has been addressed at least for 
laminar flows (e.g., Simon 1955). In the case of the strongly super-Alfvenic turbulence in the ICM an additional effect comes 
due to the bending of the field lines. Field lines are bended on scales < l m f P and this affects the ion diffusion process and 
thus viscosity. An additional suppression of the viscosity might come from the effect of plasma instabilities which affect the 
ion-ion mean free path, but that are not considered in this paper. Then, as we are interested in the compressible motions 
we discussed their generation in super-Alfvenic and MHD turbulence along with providing the estimates for collisional and 
collisionless damping of such motions. The outcome of § 3 is a validation of a basic features of scenario according to which 
the energy can be injected due to cluster mergers on large scales and energize the particles in the ICM. 

The quantitative treatment of the particle re-acceleration requires a much more rigorous treatment of mode spectrum 
and damping which was non trivial. In § 4 we make use of collisionless physics and quasi-linear theory and derive general 
formulae for the spectrum of the compressible modes, basically the ratio between the energy in magnetic fluctuations and 
the total energy in the mode, and for the damping rate in magnetized plasma, and re-obtained expressions known in the 
literature, as a particular cases of our approach. The importance of the derived formulae goes beyond our particular case of 
study, as a rigorous description of damping is important for many other astrophysical important situations, e.g. in galactic 
environments (see Yan & Lazarian 2004). 

Having at hand a description of compressible super-Alfvenic and MHD turbulence with specified injection and damping 
scales we studied in § 5 proton and electron stochastic re-acceleration by compressible modes. We focus on particle acceleration 
from magnetosonic modes and neglect the contribution from slow modes and Alfven modes. Slow modes are sub-dominant 
for particle re-acceleration as they have a phase velocity << than that of fast modes and sound waves, in addition both 
slow modes and Alfven modes get anisotropic at small scales (if injected at large scales) and this reduces the efficiency of 
gyro-resonance acceleration. We showed that because of efficient damping of fast modes at small scales the acceleration by 
gyroresonance is suppressed, i.e. only extremely high energy protons with large gyroradius can find magnetic perturbations to 
resonate with. Thus we study stochastic re-acceleration by both non-resonant large scale compressions and resonant TTD, 
and clarified the regimes when the non-resonant large scale compressions is important. The acceleration picture that is drawn 
from this paper is complex. In the case of super-Alfvenic turbulence in the ICM the turbulent bending of magnetic field 
lines limits particle spatial diffusion. Because line bending is associated with turbulent compression fast particles diffusing 
through the compressible turbulent eddies may experience efficient stochastic acceleration via Fermi 77 nonresonant -turbulent 
compression. The same particles can also experience coupling with these compressible eddies via TTD resonance which is 
found to be efficient in the ICM provided that particle pitch-angle isotropization is maintained. 

Finally, in § 7, we apply our results to the case of the particle re-acceleration scenario which is proposed to explain radio 
halos (and hard X-ray tails) in galaxy clusters. Our calculations showed that the acceleration of energetic particles in galaxy 
clusters may be efficient. Relativistic electrons in the ICM can be re-accelerated against radiative and Coulomb losses up to 
energies of several GeV (or more) assuming that compressible turbulence at large scales stores a non-negligible fraction of the 
thermal energy, namely (Vl/c s ) 2 > 0.15. These electrons would emit Mpc-scale synchrotron radiation up to GHz frequencies 
(or more) provided that the magnetic field in the ICM is at w /j,G level on these scales. In addition it also comes out that the 
re-acceleration of these electrons happens without transferring too much energy to protons, which might alleviate possible 
problems of earlier re-acceleration models that appealed to Alfven modes. 



7.2 Simplifying assumptions 

In other words, the proposed re-acceleration scenario, which makes use of compressible modes, is a plausible one and deserves 
further studies. At the moment it includes several simplifications. In particular, plasma instabilities can decrease further mean 
free path of protons, which would decrease damping of turbulence. As a result, compressible modes could cascade to smaller 
scales, making, for instance, gyroresonance acceleration by fast modes more efficient. An effect related to the gyrokinetic 
instability in accelerated particles (Lazarian & Beresnyak 2006) may act in a different direction suppressing compressible 
motions at small scales, however on the other hand the Alfvenic component generated by this instability may also accelerate 
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particles. Plasma instabilities might also affect the diffusion of fast particles and this might be important in the calculation 
of the efficiency of the acceleration from nonresonant compression. 

In addition, reconnection processes taking place in the magnetized plasma should be able to accelerate particles on their own. 
Within small volume current sheets, the percentage of accelerated particles is small. However, stochastic reconnection model 
in Lazarian & Vishniac (1999) allows acceleration of a substantial part of particles at the expense of the magnetic energy in 
the turbulent plasmas (Gouveia Dal Pino & Lazarian 2005). Therefore we believe that our treatment would underestimate 
the actual acceleration; further research should clarify the actual picture. 

Our derivations of the damping rates are valid when the ratio of the imaginary to the real part of the mode-frequency is 
much less than unity. This is generally true in the ICM and is a natural assumption for dealing with turbulence cascade, as 
in the opposite regime, no cascading is possible and the energy dissipates at the injection scale. 

7.3 Relation to Earlier Works 

Stochastic particle acceleration in galaxy clusters has been addressed by several papers (e.g., Schlickeiser et al. 1987; Brunetti 
et al. 2001; Petrosian 2001; Fujita et al. 2003). This work appeals to compressible motions to re-accelerate particles in the 
ICM. Earlier detailed time-dependent calculations of the problem of re-acceleration was addressed in Brunetti et al. (2004) 
and Brunetti & Blasi (2005), where Alfven modes were used for the purpose. Such an approach is adequate if, for instance, 
Alfven modes are injected by some mechanism at small scales. One possibility is that this might happen in the gyrokinetic 
instability scenario (Lazarian & Beresnyak 2006), and this provides an interesting possibility that we consider elsewhere. If, 
however, Alfven modes are injected at large scales the Alfvenic component at the scale of energetic particle gyroradius gets very 
anisotropic and interacts very inefficiently with the particles (Chandran 2000; Yan & Lazarian 2002). Thus fast compressible 
modes should be considered. These modes are isotropic and may scatter particles efficiently as we have demonstrated above. 

In some aspects the scenario suggested in the present paper for the ICM is similar to that adopted to calculate the 
scattering of galactic cosmic rays in Yan & Lazarian (2004), and to that adopted in the ICM by Cassano & Brunetti (2005). 
The present paper is a theoretical extension of the work of Cassano & Brunetti (2005) where a more simplified treatment 
of the resonant TTD re-acceleration of electrons by fast modes was used to derive the statistical properties of non-thermal 
emission in galaxy clusters. 

Yan & Lazarian (2004) discuss cosmic rays propagation in Milky Way thus focusing on MHD turbulence at scales I < I a and 
mostly low beta plasma. On the other hand, here we concentrated on the acceleration by motions at scales larger than I a and 
high beta plasma, the conditions which are relevant to the clusters of galaxies. This made our calculations of the particle-mode 
damping rates and of the particle-diffusion coefficients different from those in Yan & Lazarian (2004). In particular, we had 
to re-adopt many of the plasma results for high beta plasma and to treat differently magnetic field wandering. 

Finally, an additional new of this paper is that we not only considered acceleration of electrons by the TTD resonance, 
but acceleration of fast protons and electrons subjected to both the TTD resonance and the large scale compressions. 

8 SHORT SUMMARY 

The paper above explains the non-thermal emission observed in galaxy clusters as a consequence of electron re-acceleration 
by compressible turbulence. In this scenario turbulence is injected at the scale of galaxy mergers and cascades to small 
scales where the bulk of energetic particle acceleration happens. The turbulence is described by using the recent advances in 
understanding of MHD turbulence. The paper incorporates: 

I. A model of compressible turbulence in galaxy clusters. In this model the energy is injected at the scale of galaxy mergers 
and cascades to small scales where the bulk of energetic particle acceleration happens. 

II. Calculations of the plasma damping and energy of the mode for an arbitrary angle of wave propagation to magnetic field 
and a rather general model of plasma. 

III. Calculations of acceleration of protons and electrons by compressible motions in ICM plasma and a detailed application 
to the particle re-acceleration scenario to explain radio halos and possibly hard X-ray tails. 

Our results show that electrons obtain a substantial part of the energy transfered to the energetic particles, which fits 
well to the existing observational constraints. 
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APPENDIX A: DIELECTRIC TENSOR IN THE LONG WAVELENGTH LIMIT 



A relevant case for many astrophysical situations is that of long-wavelength modes for which it is \z a \ ~ kp/m a Q" << 1. 
This is also the case of the turbulent modes in the ICM of interest in the present paper. The dielectric tensor (Eq.16) is in 
the form : 



(Al) 



Barnes & Scargle (1973) calculated the tensor Rij from Eqs. 16-17 and under the conditions z a = k±p±/m a £lo << 1 and 
X" = (fc||P||/m Q - w7)/n° << 1. In this case, by taking into account that J„(z a ) = {z a /2) n ^ m (-z^/4) m /m!r(n + m+ 1), 
one finds (BS73) : 



Rij ~ — 27r^OJp jCt 



dpi 



where 
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2ft? 
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(a;-fc||t)||)- n +k\\v± 
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2X? 



EL 

p± J xf 



V 



- -! / £|| 



dp± 



+ 0(z 2 a ) 



+ m, 



aP\\ (v± 



df a (p) df a (p) 



■'-'II ■ 



dpi 



-)*a^-s}(A2) 



(A3) 



/ 



and where the term O(Za) comes from the contribution from the n > 2 resonances in Eq.(16). 

In an isotropic plasma, Eqs.(A2-A3) can be further simplified by introducing the total energy of species a 



£ a = 2nN a J J dp±p±dp\\f a (p)m a c 7, 
and the pressure of species a : 

P a = PZ= ivN a J J dp ± p ± d Pll f a (p)p ± v ± = if = 2ttA 



dp±p±dp ll f a (p)p ll v ll 



(A4) 



(A5) 



which is P a — N a kBT for a Maxwellian distribution of a particles. 

By integrating Eqs. A2-A3, introducing Eqs.(A4-A5), and requiring no net charge in the plasma (i.e., ~^2 a N a e a — 0), 
the components of the tensor are given by : 



47rfcjjc 2 

R12 = R21 = 

R13 = R 3 i=0 



UJ 
fcllC 



^{U + P a ) 



8TTk 2 ± c 2 



irk 2 ± N a 



R,, = R U ,-L^^ Pa ^^2^ {p±v i )c 



R23 = —R32 = — 



4n ik±c Na( , aL0 J J p± d pj 



R33 



^ Nael J J 



p±dp±dp\\ 



27T 2 fc||fcic i 

m 

p±v± df a {p) 
u> - fc||U|| dp\\ 

dUp) 



p±dp±dp\\ 



Pjjrj 9f a {p) 
w-fe||V|| 9pn 



dp,, 



(A6) 

(A7) 
(A8) 

(A9) 
(A10) 
(AH) 



which give all the components of the dielectric tensor when inserted in Eq.(Al). 

As usual, the integrals in Eqs.(A9-All) can be calculated using the Sokhotskii-Plamelj formula by taking into account the 
causal condition (e.g., Melrose 1968) : 

1 „ 1 



- ' JtS^LO — fc||U||^ 



w — k\\v\\ + iO uj — k\\v\\ 
where V is the Cauchy principal value, and iO is an infinitesimal imaginary term. 



(A12) 



APPENDIX B: ENERGY OF THE MODE 

The energy of the mode in a collisionless plasma is given by (e.g., Melrose 1968; Barnes 1968): 
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Figure Bl. The function F is reported, as a function of the mode-propagation angle 9, for electrons (left panel) and for protons (right 
panel). Results are shown for different values of c^/v 2 A (= f3 p i/2): 100 (solid lines), 3 (dotted lines), 1 (dashed lines), 0.1 (long-dashed 
lines). In the calculations kgT = 8.6 keV is assumed. 



B ki B ki + E k * — (uKi^j E k 



(Bl) 



where stands for the Hermitian part of the dielectric tensor (Al). In the weak damping limit (i.e. Im(uj) << 1) the 
Hermitian part of the dielectric tensor (Eq. Al) can be expressed as : 



(B2) 



where the components of the tensor Mi 



(Rfj + R°' i )/2u> 2 are given by Eqs.(A6-All). We note that the x-component of 



the electric field associated with the mode (and its spatial Fourier transform) is =0 (Eq.13), and thus only the components 
K23, K32, and contribute to Eq.(Bl). The Hermitian part of the relevant components of the tensor Mij are obtained 
from Eqs.(A9-All, & A12) and Eq.(B2) : 



TvN a (p ± v 3 ± ) 



M%3, a = -M$ 2 , a 



fcn 



. 2 k±N a ce a . 

-47T — % 

BoiOr 

dp±p± V 



dp±p_i_v± V 



dp\\ df(p) \ 

LO r — fc||U|| dp\\ J 



\ ^ 2 fc|| kjN a c 2 r 
d P\\ df(p) \ 

U!r~k m dp\\ ) ' 



dp±p±v± I V 



,(B3) 



9 P\\ J, 



The energy spectrum of the mode (Eq.Bl) is thus given by 



B k \ 2 + \E k \ 2 - 



^i a |j(«E^.«.)-^ii a |j(«E^.«.) 



Ej I £2j\\ — 1 1 



(B4) 
(B5) 

)}(B6) 



Here it is necessary to evaluate the ratio between the parallel and perpendicular fluctuations of the electric field in Eq.(B6). 
In the low amplitude regime the Fourier-Laplace transform of the mode electric field satisfies \E> • E(k,w) = 0, where ^ij = 
Rij + c 2 kikj — (k 2 c 2 — u) 2 )8ij is the Maxwell operator (e.g., Melrose 1968; BS73; Schlickeiser 2002). Thus since it is E kl — 0, 
one has : 

g|| #12 + #22 + #32 

E± - 



(B7) 



#13 + #23 + #33 

A dimensional analysis of Maxwell operator in case of the long-wavelength modes gets #32 >> #22 >> #12 and #33 >> 
#23 >> #13 (BS73), and thus the ratio between the perpendicular and parallel component of the mode electric field is : 
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Figure B2. The function V(. ■ -)/fc|| is reported, as a function of the mode-propagation angle 0, for electrons (left panel) and protons 
(right panel). Assumptions and line-styles are the same as in Fig.Bl. 

E ± ~-E^ (B8) 

which can be used in Eq.(B6) in combination with Eqs.(AlO) and (All). Here it should be noticed that it is |i?x| 2 /|-£'||| 2 ~ 
\Rs3\ 2 /\Rs2\ 2 >> 1 and thus that the fluctuations of the electric field are perpendicular to B , and to the fluctuations of 
the magnetic field (see Eqs. 12-14). On the other hand, this does not immediately imply that the contribution from |-E±| 2 in 
Eq.(B6) dominate on that from |-E||| 2 : a dimensional analysis of the elements of the tensor Rij indeed shows that |i?x| 2 A422 
is of the same order of \E\\\ M-33. 

The most important contribution of particles to the spectrum of the mode in the ICM (as in many other astrophysical 
cases) is provided by thermal electrons and protons which dominates the energy budget of the plasma. In this case the particle 
distribution function of electrons and protons is Maxwellian : 

U(P) - ^mp) - (27r)3/2 {makBT)3/ 2 ^> 
and, from Eqs.B3-B5, the components of the Hermitian part of the dielectric tensor become : 

_ C \ 2 /8VP a ivN a (j> ± v 3 ± ) \ A% 2 

-) vm Bj^) + — IaM22 + -^r Ia (B10) 

(Bll) 
(B12) 

where T a stands for the Cauchy principal value in Eqs.(B3-B5), and we put : 
16^ 2 k ]{ kjc 2 N a /— — — 

A *> = (2^)3/2 Bimi V m ° k * T > (B13) 
a _ 8ty 2 N a e a k±c , . 

(2 W )3/» Bom 3/2 fc V» T l/2' (BW) 

and 

8tt 2 N e 2 

A ™ = J^nfT 2 (m a k B T)3/ 2 k }l ' (B15) 
In the case of Maxwellian distributions the Cauchy principal value in Eqs.(B10-B12) reads : 



-M22.C 




•M23.C 


^23 

= iT.ct 




M h 33 , a 


= A 33 I a 
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Figure B3. Left Panel: The expressions Ai/ (bottom) and A2 (top) are reported as a function of the mode-propagation angle 8. 
Right Panel: The function C is reported, as a function of the mode-propagation angle 9, for electrons (thin lines) and for protons (thick 
lines) . 

In both panels, assumptions and line-styles are the same as in Fig.Bl. 

l a =(v r dP]]P J exp{- P [ _ }) =-V2^(m a k B T)i 1 ^\l-* a F(zA (B16) 
V J-oo LOr ~ k \\ v \\ 2m a k B T J a fe|| L J 



where we define 



2u a - ±ul + ^ujI for Cj » 1 



F(Cj) = 2exp{-£^} / dxexp{x 2 }^{ (B17) 
Jo 1 + 3 for d)<< 1 

which for real argument is F(x) — —1ZZ(x), and 

Z (x) =±r rftexp{ - f2} (bis) 

is the well known plasma dispersion function (Fried & Conte 1961; Melrose 1968; Percival & Robinson 1998). The adimensional 
frequency, u>, in Eq. B17 is defined (from Eq4) as : 



1/2 



^^OT = Vtfe)fe) (^) + TOFj ■ (m9) 

For d} a ~ 1 the bulk of thermal particles of species a undergoes n = O-resonance with the mode. The value of Cj a increases 
with increasing 6 and goes to infinity with 9 — * 7r/2. For a given (#, /3 P () the value of the adimensional frequency of electrons 
is about 40 times smaller than that of protons, thus electrons experience n = O-resonance with the mode at larger angles 
than protons. The value of the adimensional frequency also depends on the P p f. Cu a increases with decreasing /3 p i, while 



So — » y / 5/3y / m ct /m p / cos(#) for large (3 p i. 

The behavior of the Cauchy principal value (Eq. B16) and of F(Q a ) is driven by the value of Cb a . In Fig. Bl we report F(Cj) 
for electrons and protons for different values of /3 p i. F peaks at u> ~ 1 and for electrons this happens at larger 9 than for 
protons. With decreasing j3 p i the phase velocity of the mode increases and becomes significantly larger than the sound speed. 
This causes a shift of the peak of F toward smaller 9 in Fig. Bl, and also prevents the n = O-resonance of the bulk of the 
protons in the case of small /3 p i. In Fig. B2 we report the Cauchy principal value (Eq. B16) of thermal electrons and protons. 
The principal value goes to zero for 9 — > 7r/2 which essentially means that in this limit there is no particle-mode coupling via 
the n = O-resonance, and this is because formally infinite particle's velocity is requested to resonate at these angles. Also in 
this case the features of the curves are shifted at smaller angles with decreasing f3 p i . 
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Figure B4. Left Panel: The expression k^_S(..) is reported as a function of the mode— propagation angle 9, for electrons (thin lines) 
and for protons (thick lines). 

Right Panel: The expression fc^(2 — S) is reported as a function of the mode-propagation angle 9, for electrons (thin lines) and for 
protons (thick lines). 

In both panels, assumptions and line-styles are the same as in Fig.Bl. 

Given these results, the spectrum of the mode (Eq. B6) in a magnetized-collisionless plasma can thus be obtained in 
explicit form. After some tedious algebra, from Eq.(B6), Eqs.(B10-B17), and Eq.(B8) we find: 



W(k,oj) = 



16tt 



1 + 



(iv/kc 



-(*)" 



1 + 



N a m a c 2 +P a I] 



F a fk±C\ 
2/87T V U ) 



1+0J 2 { 33) 



1 - Hoc 1 - 



N a {p ± v'i) c 



where we put 



(l-^F Cp ) 2 (l + ^ 



{i-u> p F Qp y 

where A£l> 2 _ e = (D 2 — a; 2 , and 
A?. e 2 Bl 

{ Aj 3 ) = (k±cy(k° B Tyki Ml3vl ' e) ' 



2A 2 (P P i,0) J 8P a c 2 



}]} 



+ 



¥| 1 (4) 2 exp{-2^}(l-(^)Uxp{A^_ e }) 2 



and where 

A c m « 1 - u, a F aa (1 - a a ) + 2a a al{\ - u, a F a 

AqOW) c^J) 

with a a defined as : 
o e = -1 - 2Ai 
and 

o p = -1 + 2Ai 

in the case of electrons and protons, respectively, 



(B20) 
(B21) 

(B22) 

(B23) 

(B24) 
(B25) 
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d-^F^y i 



and where 



f -£ (X) exp{-2^} l-(^)exp{2A^_ e } 



(l-£ p Fb p ) 2 



+ ¥^(^)W^}(i-(^)Ux P {A^_ e })' 



C a (/3 p i,0) = l + 2u>l(l-u a F(u a )j 



(B26) 



(B27) 



The terms Ai, A2, and C a are reported in Fig. B3 for different values of /3 pi . 

From Eq.(B20) with (lu 2 (A 2 3 / 'A\ 2 )) » 1, and lu — V p hk, one gets the expression for the ratio between energy density in 
the magnetic field fluctuations and total energy density of the mode which is used in this paper : 

1 2 r a. /T/. x 2 q,u.\2/ \ 1 /T/.\2 



IB* 



w(k,e) 

where 



S((5 pl ,6;{f rel (p),T})= £*(Ai,0) (1 



2A 2 {P v i,e) 



8P a c 2 



(B28) 



(B29) 



is reported in Fig. B4 for different values of (3 p i \ the expression for \Bk\ 2 /W (Eqs. B28-B29) is important in our calculations 
since it allows to obtain the value of the TTD-damping rate and of the TTD-acceleration efficiency (§. 4.3-5.1). 
Both C and 5 goes to zero for 9 — > 7r/2 (Figs. B3 & B4) which means that there is no contribution to the energy of the 
mode via n = O-resonance. Also at small the contribution to the energy of the mode due to particle-mode coupling via 
n — O-resonance goes to zero (Fig.B4), and this is because in case of parallel propagation of the mode the compressible part 
(parallel) of the magnetic field fluctuations which drives the n = O-resonance goes to zero, B3 — ck±E±/uj (Eq. 14). Finally, 
the same arguments used to comment Figs. Bl & B2 can be used here to explain the evolution of the behavior of C and S 
with /3 p i. 



APPENDIX C: DAMPING COEFFICIENTS 

The damping coefficient of the modes can be obtained by the standard formula for the linear growth rate of the modes in the 
weak damping approximation (e.g., Melrose 1968; BS73): 



\ 16tyW J Ui=0 



(CI) 



where Kfj stands for the anti-Hermitian part of the dielectric tensor which can be directly obtained from Eqs. (16-17); here 
we closely follow the approach in BS73. From the Sokhotskii-Plamelj formula (Eq.A12), from Eq.(17) one has (e.g., BS73): 
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By making use of the properties of the (^-functions, from Eq.(Cl), Eq.(16), and Eq.(C2) one has 



16ui r W \k\\ I 



+ 



fell 



(C3) 



where we define 
=2 



Jl l (z a )E ± + ^-.J n (z a )E ll 



(C4) 



Here we focus on the TTD case, n = 0, which is the most important resonance of long wavelength (z a << 1) fast modes and 
magnetosonic waves. In this case it is : 



Za ^ 1 2 



n En ^ t ^l E± 2 



p± 



2Q 



(C5) 



and by using the properties of the 5-functions in Eq.(C3) we obtained a general formula for the damping rate (TTD) with 
particles of a-species : 
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(C6) 



pii ( res ) 



where 7i(:r) is the Heaviside step function (1 for cc > 0, and otherwise), the derivative of the particle distribution function 
should be evaluated at : 
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, 1/2 

(C7) 
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\ k u c ) 



and cr a is given by 



V fc ll c / V fc ± c /c« l-(^/(fc||c)) 2 V Pi 



(C8) 



Since the fluctuations of the electric field are essentially perpendicular to B (Appendix B), it is E\\ — > and a a — > and 
thus from the Faraday low (Eq. 14) the expression for the damping rate gets simplified : 
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Thus from Eq.(C9) one obtains formulae appropriate for the case of the ICM: 
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in the case of thermal particles, and 
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in the case of relativistic particles (see also BS73), where in obtaining Eq.(Cll) from Eq.(C9), one takes 
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which implies p± = p^/l — (uj/k^c) 2 , and the derivative of the distribution function is taken : 
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